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A survey is presented of phenomena in stratified and rotating fluids. Care is taken to define important 
quantities, to discuss basic concepts, to derive the fundamental equations and to present the basic nondimcn- 
sional parameters associated with these flows. After a rather extensive section on waves, the analogy between 
stratified and rotating flows is discussed. Then nonlinear processes and transport and diffusion processes are 
reviewed. Although this report is rather brief in parts, it displays the rich variety of phenomena in stratified 
and rotating fluids. It also tabulates many of the important references. 

This report also contains a discussion of some applications. Physical oceanography, physical limnology 
and meteorology are all areas of application in the earth sciences. The variety of application is stressed at the 
expense of depth and completeness. The more technological application of plumes in bodies of water or in the 
atmosphere is also discussed. 

Finally, some important problem areas arc briefly surveyed. These areas arc turbulence, nonlinear processes, 
numerical computation of flows, and further applications. 
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1 . Introduction 

An increased interest in geophysical phenomena and greater public concern about the 
quality of the environment have led to an increase in activity in geophysical fluid dynamics. 
Although fluid-dynamical phenomena are only one aspect of such broader geophysical and 
environmental considerations, they constitute a very fundamental aspect. The atmosphere, 
after all, is the layer of fluid surrounding the earth in which we live and in which the flow and 
transport influence every aspect of our lives [l]. 1 Also, the oceans and lakes are fluids covering 
the major portion of the surface of the earth, and these also influence our living [2]. Even the 
dynamics of the solid surface of the earth is now thought to be dominated by convective motions 
within the earth, and these motions are being described by the plate-tectonic model of the 
earth [3]. 

Fluid phenomena within the atmosphere, oceans and lakes are often profoundly affected 
by the variations in density within the fluid [4, 5, 6]. Such fluids, where density varies as a 
function of depth, are called stratified fluids. As a result of the increased interest in geophysical 
flows, the study of stratified fluids has been re-awakened in recent years [7, 8]. In a stratified 
fluid, buoyancy plays an important role, the stable configuration lor the fluid being one in 
which the lighter, less-dense fluid overlays the heavier fluid. The ocean-atmosphere interface 
is an example of a stratified flow; here the very large density difference between the air and the 
underlying water is such that the inertial effects of the air are usually neglected in comparison 
with the inertia of the water. However, very generally, temperature variations and chemical- 
composition and solid-particle variations provide density stratifications which lead to subtle 
and often profound differences in the flow and transport phenomena. 

A second effect within geophysical flows which can profoundly influence the flow and 
which has also undergone a resurgance of interest is the effect of rotation of the earth [9, 10]. 
The influence of rotation generally applies to larger scale phenomena than those associated 
with stratification, but often flows are affected by both rotational and stratification effects 
simultaneously [6, 11, 12, 13]. The effects of stratification and rotation are not only important 
in geophysical applications, but also in astrophysical applications [14, 15]. 

There is an analogy between the effects of rotation and those due to stratification [16]. 
Kotation of a homogeneous fluid about the vertical axis inhibits the flow in horizontal planes. 
Similarly, in a stable, density-stratified fluid, motion of the fluid in the vertical direction 
(parallel to the direction of gravity) is inhibited because of the buoyancy variation of the fluid. 
Under certain conditions the analogy reduces mathematically to the same problem. For wave 
phenomena the analogy has been known for some time, and a mathematical formalism has 
been developed to analyze these waves (as well as those arising in other flows such as magneto- 
hydrodynamic flows [17]). 

In this paper an exposition is given of stratified and rotating phenomena. There are several 
excellent texts which have come out recently (refs. [4, 5, 6, 7] and [9]) which examine stratified 
or rotating fluids, and this paper will rely heavily on these texts. However, additional features 
will also be considered. 

In the second section general considerations, such as definitions of terms, basic concepts, 
basic parameters and the fundamental equations, in the study of stratified and rotating fluids 



1 Figures in brackets indicate the literature references at the end of this paper. 
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will be presented. In the third section an examination of small disturbances (waves) to equi- 
librium configurations will be discussed. Because of the interest of the author in this area, this 
section will be rather extensive. In section four the analogy between stratified and rotating 
fluids is discussed in some detail. Nonlinear features of stratified or rotating flows are presented 
in section five, and transport and diffusion processes in these flows are surveyed in section six. 
Section seven deals with applications of stratified and rotating flows to oceanography, limnology, 
meteorology, and thermal discharges. Because the fields of application are extremely broad, 
this section will not be very detailed. Finally, section eight will present some important prob- 
lems in the study of phenomena in stratified and rotating flows, and section nine will give a 
summary and conclusions. 

2. General Considerations 

In this section a discussion is presented of the basic fluid-d}mamical considerations for 
stratified and/or rotating fluids. This discussion will attempt to be simple and clear and yet to 
introduce all of the important aspects of these systems. Reference to the several excellent texts 
on these subjects will be made frequently [4-7, 9]. 

2.1. Definitions and Basic Considerations 

For all systems of fluids the basic relations which determine the dynamics are the conser- 
vation of mass, momentum and energy. These relations are then supplemented by others which 
relate various fluid properties such as the equation of state, a relation between stress and strain 
rate, etc. The equations which express these relations will be presented in another subsection. 
The purpose of this subsection will be to define fundamental quantities and to introduce basic 
concepts related to stratified and rotating fluid systems. 

a. Equilibrium 

A stratified fluid is one in which the fluid is not homogeneous throughout, but in equi- 
librium resides in layers or strata, each of which has different physical properties. The physical 
propert}^ which most affects the dynamical behavior of the fluid is the density. Yih [7] points 
out that the density stratification between layers has two important effects, the inertia effect 
and the gravity effect. Since the inertia of a fluid element is proportional to its density, a density 
variation produces an inertia variation. 

When a body force (gravity) is present within the fluid, then the variation of density 
also affects the body force on different fluid elements. In equilibrium because of the downward 
gravitational body force, a fluid will establish itself in a stable configuration so that the fluid 
of greatest density is at bottom while the least dense fluid will be at the top. The gravity effect 
of the stratification will be the feature of stratified fluids which is emphasized in this paper, 
but inertia effects will also be discussed. 

In a stably stratified fluid, displacement upward of a small fluid element from its equilib- 
rium position leads to a downward gravitational restoring force because the density of the 
parcel is greater than that of the fluid at its displaced height. Therefore, the fluid parcel will 
be heavier than its surroundings and will tend to sink back toward its equilibrium height. 
Similarly, displacement downward of the parcel leads to it being less dense than its surrounding 
fluid, and an upward buoyancy force tends to restore it. Because of this restoring bouyancy 
force, the fluid is stable. 

The argument given above applies to a fluid in which the density does not change when 
displaced, an incompressible fluid. An essential concept when considering a compressible fluid 
is that of potential density (or potential temperature). When a fluid element is displaced, it 
experiences a change in pressure, and, corresponding to this change in pressure, a change in 
density (or temperature). If no dissipative mechanisms act, the change in pressure results in an 
adiabatic change in density (temperature). The potential density (temperature) is defined as 
the density (temperature) resulting when the fluid is compressed adiabatically to a standard 
pressure. For the atmosphere, a perfect gas is a good approximation to the equation of state 
and the potential density (and temperature) is important because of the compressibility of 
the gas. (See reference [4] for a more detailed discussion of potential density and temperature.) 

In a rotating fluid the equilibrium configuration is likewise simple, a state of rigid-body 
rotation. In this case the fluid may be described in a reference coordinate system which rotates 
with the fluid, and in this reference frame the fluid is at rest. The fluid in rotation is also stable; 
here the angular momentum of the fluid provides a mechanism for stabilization. When a fluid 
element is displaced outward from its equilibrium radius, its angular momentum will remain 
constant. Since it will have a smaller angular momentum than the surrounding fluid at its dis- 
placed radius, the Coriolis (or centrifugal )force will tend to restore it to its equilibrium radius. 
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Similarly, if the fluid element is displaced inward, its excess angular momentum will again 
tend to restore it to its equilibrium radius. 

b. Nonequilibrium 

When a stratified or rotating fluid is altered from its equilibrium state, a variety of non- 
equilibrium processes can occur. Although these processes can perhaps be classified in other 
ways, in this survey they will be classified as nondissipative and dissipative. Nondissipative 
flows are ones in which the effects of viscosity and thermal conductivity, for example, can be 
neglected. Turbulent mixing, which behaves as a diffusion process, will also be excluded in 
consideration of such flows. Dissipative flows, therefore, are ones where some dissipative 
mechanism determines the structure of the flow. Boundary layers are a prime example of such 
flows. 

Nondissipative flows will further be broken down into linear flows (or waves) and nonlinear 
flows. In a very interesting and penetrating review article [18], Ligh thill has surveyed the dif- 
ferent types and the nature of waves in fluids. The study of waves is very important for two 
reasons. First the study of waves in an otherwise steady state is a study of linear phenomena 
and linear analysis is nighty developed compared with nonlinear analysis. Therefore, most of 
the information we acquire mathematically about a system is from the analysis of small dis- 
turbances. Second, and not entirely unrelated to the first, is that essential plrysical mechanisms 
of a system come to light in the study of waves. Some nonlinear phenomena will also be consid- 
ered in this survey. 

In an inviscid, incompressible fluid with no body forces, waves can be sustained only in 
the presence of a preexisting flow or in a two-component system [18]. When a gravitational body 
force is present in a quiescent homogeneous fluid, weaves can exist at a surface only; classical 
gravitational surface waves are the examples of these waves. In a stratified fluid, however, 
waves can be sustained for which the amplitude is largest within the fluid and not at its surface. 
Such waves are called internal waves. Likewise, because of the effective body force (Coriolis 
force) in a rotating fluid, waves can be sustained internally to the fluid. Such waves are known 
as inertial waves. In this survey internal and inertial waves will be discussed in some detail 

2.2. Basic Parameters 

a. Parameters in a Stratified or Rotating Fluid 

In a stratified fluid the basic parameter which describes the stability of the fluid is the 
Brunt-Vaisala frequency. As discussed in the previous subsection, a fluid element displaced 
either upward or downward from equilibrium experiences a restoring force. This force is directly 
proportional to the acceleration of gravity g and to the equilibrium density gradient 



dz 



-t^> where p 



is the equilibrium density distribution (only a function of z), z is measured positive upward 
and gravity acts downward. The Brunt-Vaisala frequency N is proportional to the square root 
of this restoring force and is the frequenc}^ with which the fluid parcel would oscillate about 
its equilibrium position if initially displaced. The expression for N is 






N=\-—g (1) 

\ Po / 

When there is no density variation with depth (a homogeneous fluid), the Brunt-Vaisala 
frequency vanishes and the fluid is neutrally stable. 

In general the Brunt-Vaisala (B-V) frequency is a function of depth within the fluid. In 
the ocean, for example, the B-V frequency has a very characteristic behavior as a function of 
depth, being rather small in a well mixed upper layer and at great depths, but rising to a sharp 
peak generally at a depth of some hundreds of meters [5, 6]. The depth at which the B-V 
frequency is maximum is called the thermocline and basically divides two layers of fluid in 
each of which the fluid is nearly homogeneous. In large lakes the B-V frequency is usually 
similar in form to that in the ocean, with the thermocline occurring at much smaller depths 
[6, 19]. 
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The basic stability considerations are changed somewhat when the fluid is compressible, 
as discussed in the preceding subsection. In a compressible fluid, acoustic as well as internal 
waves can be sustained within the fluid, each type of wave depending upon different 
mechanisms. A basic parameter which arises in the dynamics of a compressible fluid is the speed 
of sound c, the speed with which an acoustic signal propagates [20]. When the compressibility 
of the fluid is also taken into account in a stratified fluid, the Brunt-Vaisala frequency must 
be modified to be 



N= 



■Ai+% 



(2) 



In the atmosphere, where the compressibility must generally be accounted for, the B-V fre- 
quency generally passes through two minima and two maxima with increasing altitude [6]. 
In a rotating fluid the corresponding parameter is the angular velocity 12 of rotation of the 
fluid. In this case a fluid parcel displaced from its equilibrium radius will oscillate at twice the 
angular velocity about its equilibrium position. Therefore, the parameter 212 serves the same 
purpose within a rotating fluid as the stability parameter N does in a stratified fluid. 

b. Dimensionless Numbers 

Several nondimensional parameters arise in the study of stratified and rotating flows. 
These parameters can be divided into those which describe local properties within the flow 
(and therefore vary from point to point within the fluid) and overall dimensionless numbers 
which characterize the whole flow. (Batchelor [21] discusses this difference in some detail.) 
In this paper mainly the latter type of dimensionless parameter will be discussed. (See refer- 
ences [4] and [9] for more complete discussions of dimensionless numbers in stratified and 
rotating fluids.) 

To define overall dimensionless numbers it is necessary to select characteristic scales of 
velocity U, length L, density p and density variation Ap for the flow. (With this notation 
an overall Brunt-Vaisala frequency can be defined as 

AT 2 _#Ap _ 



An obvious dimensionless number characterizing the stratification of the fluid is the density 
ratio n=Apo/pQ. 

A second dimensionless number, which will be used in this paper is the Richardson num- 
ber, defined as 

The Richardson number compares the buoyancy to the inertia terms in the momentum con- 
servation relations. When the motions are steady, the effect of variations in density in the inertia 
terms can be ignored (the Boussinesq approximation) and the flow is nondissipative, then the 
Richardson number is the only dimensionless number characterizing the flow [4]. 

Often the Froude number is used to characterize surface flows; the Froude number is defined 
as 

Fr= • (4) 

{gL y,2 W 

It compares a characteristic flow velocity to the velocity of long waves on a free surface [4]. 
For stratified flows the internal, or densimetric Froude number, defined as 



F -7^-™ (5) 



'gtf 



is often used instead of the Richardson number. The densimetric Froude number is the inverse 
square root of the Richardson number. In this paper the Richardson number will be used con- 
sistently. 
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When the flow is dissipative, additional dimensionless parameters can be defined because 
additional physical effects also become important. In a flow where viscosity is important, the 
Reynolds number is also a relevant parameter. 

ReJ^ (6) 

where v is the kinematic viscosity. As is well-known, the Reynolds number compares the inertia 
to viscous effects. In geophysical flows, it is generally very large, indicating that boundary 
layers may be important and that turbulence is probably a dominant feature in these flows [5]. 
Another fundamental number is the Prandtl number, defined as 

Ft=v\k (7) 

where k is the thermometric diffusivity. The Prandtl number compares the viscous diffusion of 
momentum to the diffusion of heat. 

In a rotating fluid the key parameter in the absence of dissipative effects is the Rossby 
number, defined bv 

-&■ 

It is the ratio of the convective acceleration to the Coriolis acceleration and also provides an 
estimate of the importance of nonlinear effects [9]. When e<<l, linear theories are appropriate, 
whereas nonlinear analysis must be used when e approaches unit magnitude. 

Dissipative effects of viscosity are measured by means of the Ekman number, defined as 

E =m (9) 

which compares the viscous force to the Coriolis force. It is generally quite small, indicating 
the potential importance in rotating fluids of boundary layers. 

When both rotation and stratification are important, both sets of dimensionless param- 
eters must be considered. Other nondimensional numbers can be formed from these basic ones 
and will be introduced as needed. 

When time-dependent phenomena, are considered, a frequency of oscillation co is also a 
relevant parameter. Then two additional nondimensional numbers o)/N and co/2^2 will also be 
important. The first is the ratio of the frequency of oscillation within the flow to the character- 
istic B-V frequency while the second is the ratio of the frequency of oscillation to twice the 
rotational frequency of the fluid. These ratios will be important when wave propagation is 
considered in section three. 

When the compressibility of the fluid is important, the velocity of sound c also becomes 
a relevant parameter. General^, the sound speed varies from point to point within the fluid. 
A nondimensional number describing the flow with compressibility effects is the Mach number 

M=U/c (10) 

where c is a characteristic sound speed (say the maximum value). Generally, in stratified and 
rotating fluids the compressibility is neglected. 

2.3. Typical Magnitudes of Characteristic Parameters 

In this subsection a tabulation is presented of the magnitudes of the parameters discussed 
in the preceding subsection. These estimates provide a guide to analyses of stratified or rotating 
flows since they indicate which physical effects can be expected to be important in various flows. 

Table 1 presents for reference the values of the acceleration of gravity and the angular 
velocity of the earth. Then a list of representative values of the basic physical properties for 
air and water are tabulated. Since these properties are functions of temperature and pressure, 
in general, it is necessary to determine the temperature range and pressure range of interest 
before appropriate values can be determined. These values have simply been listed for reference. 

Also listed in table 1 is a range of values of the Brunt-Vaisala frequency and period for 
the ocean, for lakes and for the atmosphere. From this tabulation it is apparent that there is 
a wide variation in the values of this parameter and, consequently, a wide variation in the 
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effect of stratification. Finally, the magnitude of the density variation across the thermocline 
in the ocean or in lakes is given in this table. 

Using values in this table and values of appropriate length and velocity scales, estimates 
of the magnitudes of the nondimensional numbers can be made. First, the ratio of density 
variation across the thermocline to absolute density is approximately n^lQ~ 3 . 

Table 1. — Representative values of parameters 



Acceleration of Gravity and 


Angular Velocity of the Earth 


Parameter 


Value 


9 

2Q 
T a =2ir/Q 


981 cm/s 
1.45X10" 4 rad/s 
24 h 



Basic Physical Properties of the Fluid 



Water 




Air 


Parameter 


Value 


Parameter 


Value 


po (4°C) 

ix (20° C) 

v — fi/p 

Co 

k = thermal conductivity 

c P 

K = k/pC D 

Pr=v/K 


1. g/cm 3 

1. x 10-2 g/cm-s 

1. x 10-2 cm 2 /s 

1. 5 x 10 5 cm/s 

5. 8 x 10- 3 W/cm-K 

4. 18J/g-K 

1. 4 x 10-3 cm 2 /s 

7.2 


po (0°C, 1 atm) 

n (20°C) 
v = fx/ Po 
Co 

k 

K 

Pr 


1. 29 x 10-3 g/cm 3 
1. 81 x 10" 4 g/cm-s 

1. 36 x 10- 1 cm 2 /s 
3. 5 x 10* cm/s 

2. 4 x 10" 4 W/cm-K 
1. OJ/g-K 

1. 80 x 10" 1 cm 2 /s 
7.6 



Range of Values of the Brunt-Vaisala Frequency and Period 



N 

T N = 2ir/N 

Ap (across thermocline) 



Ocean 



10- 2 -10- 4 rad/s 
10 min-2 h 
10" 3 g/cm 3 



Lake 



5x 10-2-10" 
2 min-2 h 
10 -3 g/cm 3 



rad/s 



Atmosphere 



3 x 10- 

3 min- 



2 -10-2 rad/s 
-10 min 



Values of the Richardson number (or the internal Froude number) depend upon the values 
of the velocity scale U and the length scale L. In meteorological applications, velocities are 
typically measured in meters per second or greater and lengths are measured in meters or greater 
(kilometers to thousands of kilometers). In oceanographic applications velocities are typically 
much smaller, being measured in centimeters per second with lengths comparable to those in 
atmospheric considerations. Therefore, the Richardson number is typically rather small, rang- 
ing from values at most of order unity down to values of order 10~ 4 or smaller. 

Here, it should be noted that the vertical length and velocity scales are characteristically 
smaller than the corresponding horizontal scales. Because the oceans, lakes and the atmosphere 
are relatively shallow, but laterally extensive bodies of fluid and because the stratification in- 
hibits only the vertical motion, the characteristic vertical length and velocity are often much 
smaller than the corresponding horizontal quantities. In the Richardson number, many times 
the relevant length scale is the vertical one while the velocity scale is the horizontal one. 

In the laboratory, typical length and velocity scales are much smaller while the density 
variation is generally comparable to that in geophysical fluids [22, 23]. Because the velocity 
scale appears quadratically in the denominator while the length appears only linearly in the 
numerator, the Richardson number is characteristically much larger. If the velocity is measured 
in centimeters/second and the length in centimeters, then the Richardson number is of order 
unity. 

The Reynolds number varies even more dramatically between geophysical and laboratory 
flows. Since the kinematic viscosity varies by over an order of magnitude between water and air, 
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characteristic Eeynolds numbers will also vary between oceans and lakes on the one hand and 
the atmosphere on the other. In oceans or lakes, taking a length scale of one meter and a velocity 
of one meter per second, the Eeynolds number is of order 10 6 . In the atmosphere, with the same 
length and velocity scale, ^10 5 Re. In the laboratory, taking a length of one centimeter and a 
velocity of one centimeter per second, in water, Re=10 2 whereas in air Re^lO. 

The Prandtl number is a characteristic of the fluid and therefore has been listed in table 1. 

In a rotating fluid the Rossby number can be calculated for geophysical fluids and labora- 
tory fluids. In geophysical applications, a velocity of a meter per second and a length of a meter 
yield a Rossby number, which is very large; hence these scales are not consistent since the 
Rossby number must be of order unity or smaller [9]. The velocity must be smaller or the length 
longer to make the Rossby number of order unity; values of £7=10 cm/s and L=10 5 cm would 
be much more appropriate. In the laboratory representative values given by Greenspan [9], 
are 12=1 rad/s, L=20 cm and e=0.01 (or £7=0.2 cm/s). 

The Ekman number is of order 10 ~ 8 for oceans or lakes when a length of 10 5 cm and the 
viscosity of water is used. Using the same length scale and the viscosity of air, in the atmosphere 
the Ekman number is of order 10~ 7 . (It should be noted that many analyses are conducted 
using a much larger kinematic viscosity to represent the fact that the motion is turbulent. 
This would change the value of the Ekman number considerably, although it would probably 
still remain small.) For the laboratory, using 0=1 rad/s and i=20 cm and water, the Ekman 
number is of order 2 X 10~ 5 . Therefore, in all cases the Ekman number is very small, implying 
that viscous effects in the main body of the flow are small and that bounda^ layers will arise. 

Table 2 presents a list of nondimensional numbers calculated in this subsection. The Mach 
number, based upon a velocity of one meter per second and the sound speeds given in table 1, 
is seen to be very small. 





Table 


2. — Representative values of dimensionless numbers 


Dimensionless Number 


Value 






Lakes & Oceans 


Atmosphere 


Laboratory 


n 


10" 3 




10- 3 


Ri 




1-10- 4 




lO-lO" 1 


Re 




10 6 


10 5 


10 2 


Pr 




7. 2 


7.6 


7.2 


€ 




1 


1 


0.01 


E 




10" 8 


io- 7 


2X10- 5 


M 




<io- 2 


<10~ 2 


<10~ 2 



2.4. Fundamental Equations 

The fundamental equations are expressions of conservation laws for mass, momentum and 
energy. These equations are derived in various forms in each of the primary references [4, 5, 
6, 7] and [9]. Therefore, in this subsection, the equations will simply be stated in suitable form 
and some observations will be made concerning them. 



a. The Equations of Conservation of Mass and Momentum 

The conservation of mass for a compressible fluid is written in standard form as 



dt 



+V-(pu)=0 



(11) 



where p is the density of the fluid, u the velocity, t time and the coordinates are x h x 2 , and x 3 
(or sometimes x, y, and z will be used). This equation can be written in several alternate forms, 
and these forms will be used as needed. 

The equation for the conservation of momentum is somewhat different than the standard 
form written in most introductory texts on fluid mechanics. Because the fluid is stratified and 
gravity g acts as a body force per unit mass, the equations must include an expression for this 
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body force. In addition, when a rotating fluid is considered, this equation is written with respect 
to a reference frame that rotates with the fluid at a constant angular velocity 12. This equation is 

p^+2pOXu+pGX(OXr)+V#-pg=f. (12) 

Here, p is pressure, r= (x h x 2 , x 3 ) is the radius vector to a fluid element and f is the resultant of 
all other forces acting on a unit volume of fluid. Both the angular velocit} r 12 and the gravitational 
body force g have been made into vector quantities since each has a direction. The time deriva- 
tive is the substantial derivative 

D d , 

The two terms involving £2 in the equation represent the effective body force per unit 
volume due to the rotation. The first term is the Coriolis force while the second is usually 
called the centrifugal force. Often the centrifugal-force term is absorbed in the pressure term 
by defining a reduced pressure. (See ref. [9] for details.) 

The terms on the right-hand side will generally be taken to be a result of viscous forces 
alone. When viscous stresses are considered, it will also be assumed that the fluid is incompres- 
sible and f=p ^V 2 u will be used. Otherwise it will be assumed that the right-hand size is zero. 

b. Equilibrium 

In equilibrium the simplest case to consider is when the fluid is in a steady state so that 
^-=0, and there is no velocity of the fluid in the rotating reference frame. Then 

PofiX(S2Xr)+V^o-Pog=0. (12a) 

Two basic equilibrium states will be considered, the state of no rotation in a stratified fluid 
and the state of no stratification in a rotating fluid. 

When there is no rotation, the equilibrium is determined by this equation with the first 
term set to zero, the so-called hydrostatic equilibrium state: 

Vj)o=gPo. (12b) 

In this case the pressure gradient is parallel to the gravitational body force and then the density 
and pressure depend upon one coordinate parallel to g. 

The other basic equilibrium state is one in which the fluid is unstratified but rotating. In 
this state the last term in eq (13) is zero, and the pressure gradient is balanced by the centrifugal 
force. A great deal can be learned about the character of flows which deviate somewhat from 
this state of equilibrium by a fundamental theorem for rotating flows called the Taylor- 
Proudman Theorem [9, 72]. The Taylor-Proudman theorem states that when the flow in a 
rotating reference frame is steady, inviscid and slow (so that the Rossby number is small, 
e<<0)j then the motion does not vary in the direction of the rotation vector. Howard [72] 
has given a very clear and concise discussion of this theorem together with a discussion of other 
important ideas such as Taylor columns, geostrophic flow and deviations from geostrophic 
flow. The following presentation follows his rather closely. 

For an incompressible fluid, the equation of continuity (11) reduces to the statements that 
V-u— and p is constant following a fluid element. For slow, steady and inviscid flow, eq (12) 
reduces to the equation [9, 72] 

2QXu=-V(p/p-l/2(fiXr)-(flXr)+{/2). (13) 

The Taylor-Proudman theorem follows directly by taking the curl of this relation and noting 
that £2 is constant and v-u=0: 

fi-VU=0. (13a) 

Therefore there is no variation in the velocity in the direction of the rotation. 
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Equation (13), representing a balance between the (reduced) pressure force and the Coriolis 
force, is known as the "geostrophic approximation" [1, 6, 9, 72]. When the rotation and buoyancy 
force are taken as vertical, then the vertical component of eq (13) is simply the hydrostatic 
equilibrium equation, eq (12b). The horizontal components of this equation then yield the 
balance between pressure and Coriolis forces. 

ox 

(14) 

2V P v=- d -P 
dy 

where u and v are the x and y or horizontal components of the velocity, and p is the reduced 
pressure [9]. 

Equations (14) show that [72] 

v d /+u d /=0. (15) 

ox ay 7 

Hence the velocity vector is perpendicular to the pressure gradient in the horizontal plane, 
or the stream lines coincide with constant-pressure contours (isobars). Many large-scale geo- 
physical flows are approximately geostrophic [6]. 

If the fluid is bounded by a horizontal surface (perpendicular to the axis of rotation) as 
would be the case in the laboratory for example, then the vertical velocity must be zero at this 
surface. By the Taylor-Proudman theorem, then, the vertical velocity is zero throughout so 
that the flow is purely two-dimensional and horizontal. 

Taylor [73] demonstrated the validity of these considerations with an experiment which 
is well known now. In this experiment he slowly dragged across the bottom of a rotating tank a 
small circular bump. In the absence of motion of the bump, the fluid in the tank rotates as a 
solid body. But when the bump is dragged slowly across the bottom, all motions relative to the 
tank are two dimensional so that the column of fluid directly above the bump is dragged along 
above the bump. This column (and similar columnar flow) has become known subsequently 
as a Taylor column, and a great deal of effort has been expended to understand how 
such columns form and their structure [9]. 

c. The Energy Equation and the Thermodynamic Relations 

The energy conservation equation can be written in a variety of ways, and is usually 
written in a form suitable for the particular application. A general and convenient form is 

TJre/r, (16) 

where s is the entropy per unit mass, q is the energy added per unit mass per unit time, T is 
the temperature and D/Dt is the substantial derivative. 

The dominant form of energy addition to a fluid element will be taken to be heat 
conduction. Viscous degradiation of mechanical energy also adds heat to an element, but this 
energy addition is small [6, 9] and will not be considered. Radiative energy addition in the 
atmosphere and at the surface of water bodies is important, but will not be examined here. 
When heat conduction is important, the energy addition term is 

2 =-V-(*vD. (17) 

P 
Otherwise q will be taken to be zero. 

From thermodynamics relations eq (16) can be written in other forms; one particular form 
which will be useful later is 

DT Dp 
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where C p is the constant-pressure specific heat and 

p \dT/ p= constant 

is the coefficient of thermal expansion. As noted by Greenspan [9], the second term on the left is 
generally small in liquids. If this term is neglected and k is assumed constant, then eq (18) 
reduces to 

DT 

where 

K=Jc/pC p . 

In this equation convection and thermal diffusion are important for determination of the tem- 
perature distribution. 

Different equations of state and thermodynamic relations will be considered depending 
upon the fluid being considered. The atmosphere is a mixture of gases, predominantly nitrogen 
and oxygen; but in the atmosphere an equation of state appropriate for a single perfect gas of 
average molecular weight 29 can be used. Hence 

p= P RT (20) 

where 

R=2.S7Xl0 2 J/kg-K 

can be used with the ratio of specific heats 7=1.4 (ref. [6]). The other thermodynamic relations 
are also those appropriate to a perfect gas. 

Some limitations to the use of these equations should be noted. First, the atmosphere con- 
tains water vapor, and this water vapor, through the mechanisms of evaporation and conden- 
sation, has major effects upon the" dynamical processes. Water vapor and other constituents 
of the atmosphere, most notably carbon dioxide, also play a fundamental role in the thermal bal- 
ance of the atmosphere through absorption and re-radiation. (The temperature variation in the 
atmosphere is determined in large part by the incident solar radiation and by the low-tempera- 
ture radiation from the surface of the earth.) The atmosphere contains particulate matter, and 
this material is also important. Finally, at higher altitudes the assumptions of a diatomic, neu- 
tral gas may break down since the gas dissociates or ionizes. (At even higher altitudes the as- 
sumptions that the gas is a continuum breaks down.) However, none of these considerations will 
be taken into account in this report. 

In lakes the equation of state and thermodynamic relations appropriate to fresh water can 
be used [6]. The only point to be noted with respect to these relations is that the basic co- 
efficients in various relations are not constant. Rather, water at atmospheric pressure has its 
maximum density at about 4°C, and the temperature at which the maximum density occurs 
is a function of pressure Therefore, the coefficient of thermal expansion for water vanishes at a 
certain temperature, which itself is a function of pressure. This variation has important im- 
plications for the stratification cycle of lakes [6, 22]. The equation for the variation in the co- 
efficient of thermal expansion in pure water can be taken as [6] 

a=— (|^) =15X10-%, (21) 

where 

T M =r_4.0+2X10- 8 ^ (22) 

with p given in atmospheres and T in degrees Celsius. The density variation from a reference 
density is then taken to be proportional to the variation in temperature from the reference 
temperature with the coefficient for thermal expansion being the proportionality factor. 

In the ocean, salinity plays a very important role in the thermodynamic relations. The 
salinity is the mass of dissolved solids per unit mass of sea water. Once again, the coefficients 
in these relations cannot be taken as constants: for example, the coefficient of thermal expansion 
can be expressed as 

a=9.5X10~ 6 T' M (23) 
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where 

TV=r+7.o+2xio- 9 ^ 

with p given in atmospheres and T in degrees Celsius. 

In general, the salinity of the water is not constant from position to position, and an 
equation of state for sea water should consider salinity, as well as two other basic thermo- 
dynamic variables (such as pressure and temperature), to be a basic thermodynamic variable. 
If this procedure is followed, an equation of continuity for the salinity and other equations 
relating the flux density of salt to the temperature and salinity gradients must be added to 
the basic equations. 

When salinity gradients are considered, it can often be assumed that the density can be 
expressed as a linear function of the variation in salinity from some reference value. Therefore, 
with p ro the reference value of density and T m and S m the reference values of temperature and 
salinity, the complete variation of density can be expressed as 

p= p „[l-a(T-T n )+p(S-S 9 )] (24) 

where 

aand H(l§L 

are "the coefficients of expansion. " The energy equation will be taken to be eq. (19), and the 
equation for salinity conservation will likewise have this form 

P^-V-(pt>VS)=0 (25) 

where D is the diffusion coefficient for salinity [4, 5]. 

A brief observation should be made regarding the Boussinesq approximation. If the 
compressibility of the fluid can be neglected (as is generally the case in water and often the 
case in the atmosphere), the Boussinesq approximation can be made [5]. In this approximation, 
variations in density from a constant reference density are considered to be significant only in 
buoyancy terms (that is, terms multiplied by the acceleration of gravity) . In inertia terms the 
actual density is replaced by the reference value. This approximation is generally excellent 
and will be used except in a few parts of this report. 

3. Small-Amplitude Waves 

Ligh thill [18] has presented a survey of the different types and the nature of waves in 
fluids. In that survey he discussed the physical mechanisms for wave phenomena and some 
mathematical aspects of the formulation and solution for such waves. Furthermore, in a series 
of papers [17, 18, 24, 25] he has presented a mathematical formalism for the solution of the 
equations of a general class of these wave motions. A recent book by Whitham [26] presents a 
detailed and rather complete discussion of waves, the formulation of equations describing 
these waves and methodology for solution of the equations. 

In this section a discussion of waves will be presented which relies heavily on the papers 
and monograph referenced above. The emphasis of this discussion will be on waves in stratified 
and rotating fluids. In the first subsection general considerations of waves will be presented 
including the classification of these waves into hyperbolic and dispersive classes. In the second 
subsection the general methods will be discussed for analysis of equations describing small- 
amplitude waves. In the final two subsections waves in stratified and rotating fluids respec- 
tively will be surveyed. 

3.1. General Considerations 

Whitham [26] defines a wave as "any recognizable signal that is transferred from one part 
of the medium to another with a recognizable velocity of propagation." As general as this 
definition seems to be, it precludes light waves, which do not require a medium in which to 
propagate. Therefore, the definition used here will modify this definition slightly to read: 
a wave is any recognizable signal that is transferred from one position to another with 
a recognizable velocity of propagation. 
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Whitham then proceeds to define two main classes of waves, which he refers to as hyper- 
bolic and dispersive waves. The former he defines with reference to the form of equation from 
which the waves arise, and the latter he defines with reference to the type of solution, in the 
linear case, which the system admits. He points out that these systems are not exclusive, but 
there is some overlap between the two. 

Within each class he first discusses the linear equation from which the class gets its name. 
For the first class, the rather familiar ideas of a lryperbolic equation and the corresponding 
characteristics for this equation arise. The ideas of a dispersion relation, phase velocity and 
group velocity arise from a discussion of linear dispersive waves. Then nonlinear processes 
within eacli class are examined. For lryperbolic systems this leads to the ideas of overlapping 
characteristics, discontinuous solutions and, in the case of gasdynamics, shock waves. For 
dispersive systems the generalization to nonlinear cases leads to dispersion relations where 
the amplitude enters the relations. The analysis of nonlinear, dispersive waves is an area of 
very active research, and Whitham has been an important contributor to this area. In this 
subsection some of the basic ideas mentioned above will be discussed for linear systems; the 
present report cannot examine nonlinear waves. 

a. Hyperbolic Waves 

Waves for which the mathematical description satisfies hyperbolic partial differential 
equations are called hyperbolic waves [26]. To discuss hyperbolic equations, a few definitions 
must be given, and general considerations related to hyperbolic equations will be discussed. 
Then the specialization of these general considerations to the prototype hyperbolic equation, 
the wave equation, will be presented. The presentation of this section will follow that of refer- 
ence [26], chapter 5. 

A quasilinear system of first order equations is a system linear in the first derivatives of the 
dependent variables. The coefficients in the quasilinear system may depend both upon the 
independent variables and upon the dependent variables (making the system nonlinear). 
These systems can be examined for n dependent variables and also an arbitrary number of 
independent variables. However, the most useful system to analyze is one consisting of n 
dependent and two independent variables, taken to be time t and position x. For this system 
the general equations are 

Aij jt +a,jd S +bi=(} i=i > 2 >-..,« (i) 

where the n h j=l, 2, . . ., n are the dependent variables and A ij} a ?? and h t will generally 
depend upon u h x, and /. 

Characteristics for this system of equations are curves in the x, £-plane along which the 
directional derivatives of the dependent variables are related: that is, along a characteristic 
curve the partial differential equations (P.D.E/s) reduce to ordinary differential equations. 
Linear combinations of the P.D.E.'s can be found 



i ( A ^+^W +l ^° 



(2) 



such that all derivatives are in the same (characteristic) direction. Then the equations can 
be written 

^w+«©+^= - (3) 

The condition for such l, to exist is that 

A ti -^a,j=0 (4) 

a 

for some ratio j3/a. The curves 

t=T( v ), x=X( v ) 

for which 

dT dX ,., 

d^ = P' w =a (5) 
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are the characteristics. Here, 77 is a parameter which varies along the curve. The equation 

is then said to be written in characteristic form. It is an ordinary differential equation relating 
variations of the dependent variables along the characteristic direction. 

A definition of a hyperbolic P.D.E. can now be made. A system of equations (1), satisfying 
the condition that the matrices in the equations not be singular 

IXAij+fxa^O (7) 

for some X and /x, is hyperbolic if n linearly independent real vectors 

P\A:=i,2, . ,.,w 

can be found such that 

^*>(A«a»-a </ j8 (W )=0 (8) 

for each k, and the corresponding directions (a ik \(3 {k) ) are real with 

(a w ) 2 +(^ w )V0. 

For a definition of hyperbolic systems of quasi-linear P.D.E. extended to m independent vari- 
ables, see reference [26]. 

As an illustrative example of a system of quasi-linear hyperbolic equations, consider the 
equations of gasdynamics for one-dimensional, unsteady, isen tropic flow: 

dp dp . du 

m +u Fx +p d-x =0 



dt dx 



Here p is the gas density, u its velocity and p its pressure, c 2 is the speed of sound in 
the gas, defined by the relation 



HP- 



\dp) 



the derivative taken at constant entropy. The first equation is that of conservation of mass, 
the second is conservation of momentum and the third is conservation of entropy. The last 
equation is already in characteristic form. The first equation multiplied by c 2 , the second by 
±cp and the two added to the third equation produce the other two characteristic forms. 
Therefore, 

&-<—> i-"K +< «-°>tO=° (io) 

or 

along the characteristic defined by -r;=u J rc i j, + cp ~T f = ^ 
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along the characteristic defined by ~ff ==u ~ c J ~~h~ c P ~J+ = ® (H) 

and 

along the characteristic defined by -r;=u, -r-.—c 1 -r=0. 

* at at at 

Hence, this set of equations is hyperbolic. 

The equations of gasdynamics have formed the basis for much of the development of the 
theory of nonlinear hyperbolic equations. As can be seen, the characteristics are dependent 
upon the solution itself, and therefore the equations are nonlinear. The nonlinearity can lead 
the characteristics to overlap, and this "breaking" of compression waves produces shock waves. 
As pointed out already, the nonlinear aspects of the weaves will not be considered here. 

These equations, when linearized, lead to the wave equation, the classical hyperbolic 
equation describing waves. To linearize these equations, take as the undisturbed state a quies- 
cent fluid ^o=0 with constant properties 



Po ,PoancU 2 =(g) ; 



Then let p, u, and p represent the (small) deviation of the respective quantities from their 
undisturbed values. The equations can be written 



dp , dp^ fdu du\ 
' \dt +Co diJ = 



dt ' + Co q- x + c oPo ( xT + c o -u l = ( 12 ) 



dp dp fdu du\ n , 1oN 

~-c 0dx -c oPo ^-c 0d -)=0 (13) 

The third equation becomes unnecessary because it simply relates density and pressure fluctu- 
ations through the sound speed. Let 



operate on eq (12), 



(d d\ 

\M~ C °d~x) 

\di +Co Fx) 



on eq (13), and add these relations. The result is the one-dimensional wave equation for the 
pressure fluctuation 

(5e-*'a?>-° (14) 

The general solution for this equation is 

p=Kx—c t)+g(x+c t) (15) 

and this solution describes small amplitude waves. The first function /represents a wave which 
propagates unchanged in form to the right while the second function g represents a wave propa- 
gating without distortion to the left. When the disturbance pressure and its time-derivative 
are known everywhere initially, these two functions are determined uniquely. 

Equation (14) can also be written in factored form as 

(!-«iX! + *i) ? =°- <16) 
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Either one of the factors represents the simplest linear hyperbolic equation, each factor giving 
rise to a wave solution. For example, 

(I+ C °J:>~ =0 < 17 > 

has as its general solution 

p=j(x—c t). 
A more general, quasi-linear form of this simple and important hyperbolic equation is [26] 



(£+ c ^>=° 



(18) 



where the tilde has been dropped and c(p) is a function of the dependent variable. 

b. Dispersive Waves 

Perhaps the more interesting and certainly the more challenging class of waves at this 
time is that of dispersive waves. For the linear case, a dispersive wave is recognized by its form 
of solution (j>, having the form 

ct>=A exp (id) (19) 

for elementary solutions. Here A is the amplitude and is the phase. In the simplest case A is 
a constant and 6 has the form 

6=k.x-ut. (20) 

Here k is the wave number and o> is the frequency and each is constant. Since the equations 
are linear, the amplitude factors out and is arbitrary, while the frequency and wave numbers are 
related by the governing equations through the so-called dispersion relation 

0(«,k)=O. (21) 

The dispersion relation can be solved for the frequenc}^ as a function of wave number. In 
general there will be more than one function u=W(k) which satisfies the dispersion relation; 
each such solution is called a mode of the system. Since the governing equations are linear, 
various modes of the system can be superposed. 

When 9 is taken to be constant, surfaces of constant phase are determined. In the simplest 
case, surfaces of constant phase are planes. The phase determines the position in the c}^cle 
between a crest, where the real part of <f> is maximum, and a trough, where real <f> is minimum. 
The normal distance between two surfaces of constant phase representing successive crests 
(or troughs) is the wave length X and is equal to 

X=2ir/|k|. (22) 

The planes of constant phase are seen to change position, but not orientation as time 
increases; therefore, the surfaces of constant phase propagate as time evolves. At a fixed position 
the solution is seen to undulate as successive crests and troughs pass the point. The period r 
between successive crests is equal to 

t=— • (23) 

CO 

The propagation velocity for the surfaces of constant phase is the phase velocity c for the wave 
and is equal to 

c =±ft"t (24) 



where 



is a unit vector in the k direction. 



T 



Ikl 



fc ~|kl 
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The meaning of the term dispersion can now be seen. A wave of particular wave number k 
propagates at a specified phase velocity c determined by the frequency u=W(k) and the wave 
number. For another wave of different wave number, the freuqnency and the phase velocity 
will, in general, be different. Hence the waves will disperse as time evolves. Witham [26] defines 
a dispersive wave to be one having the form (19) with phase (20) satisfying the dispersion 
relation (21) such that, if w=W(k) is a mode, W(k) is real and the determinant 



d 2 W 
dk l dk j 



^0. (25) 



The last two conditions guarantee that the solution will in fact be wave-like and that the phase 
velocity will not be constant for all wave numbers (so that the waves do disperse). 

A very important concept in the discussion of dispersive waves is the concept of group 
velocity. The group velocity is the velocity with which a "group" of waves, or a general wave 
composed of a distribution of elementary waves of different wave number, would propagate. 
This concept can be introduced in a variety of ways [24, 26]. Here it will simply be noted that, 
in the one-dimensional case, when two waves of equal amplitude and wave numbers k { and k 2 
and frequencies a>i and o> 2 respectively are superposed, their sum can be written 

cj)=A cos (kix— uit)+A cos (k 2 x—co 2 t) 

(26) 

=2A cos (^ x-^ t) cos (*!=*' x-^ t} 

In this latter form the superposition is seen to be the product of a wave with average wave 
number (ki-\-k 2 )/2 and average frequency (coi + co 2 )/2 and a relatively slowly varying, modulating 
wave with wave number (ki — k 2 )/2 and frequency (co 2 — o?i)/2. A plot of this function for a 
fixed time shows that the amplitude of the high frequency or short wavelength wave (the first 
cosine term) is modulated by the low frequency wave over a wavelength 47r/(Z; 2 — k\). The low 
frequency or modulating wave propagates with a velocity 

which in the limit of k 2 —*ki and w 2 — »«i, becomes the group velocity in the one-dimensional case 

«<-ff <»> 

In the general three-dimensional case, the group velocity is defined by the equation 

c^=Vrco (29) 

where Vk means the gradient with respect to the components of the wave number vector k. 
It can be shown that energy and several other features propagate with this velocity [26]. 

When the medium within which a wave propagates is anisotropic, waves which are launched 
in one direction propagate with a velocity different from those launched in a different direction. 
Then both the phase velocity and group velocity depend upon the direction in which the waves 
propagate. However, in a dispersive, anisotropic medium a very surprising phenomenon arises: 
generally the group velocity is in a different direction than the phase velocity. Therefore, in 
general, energy propagates in a direction different from that in which individual component 
waves propagate, leading to surprising and interesting results. A general mathematical formula- 
tion of these waves by Fourier anatysis will be discussed in the next subsection. 

Stratified fluids, rotating fluids and electrically conducting fluids in the presence of magnetic 
fields are all examples of anisotropic media with respect to wave propagation. In each case 
there is a direction specified by the constraint on the fluid, gravity, rotation and the magnetic 
field, and waves propagating along the direction of the constraint experience different forces 
than do those propagating transverse to the constraints. The features of waves in stratified or 
rotating fluids will be discussed in some detail in remaining subsections. 
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3.2. Methodology for Analysis of Waves in Dispersive, Anisotropic Media 

As mentioned above, both Lighthill and Whitman have been very instrumental in the 
development of general methodology for the analysis of waves in dispersive and anisotropic 
media. Some of the extensions of these methods make them also applicable to inhomogeneous 
media, and considered in the broadest context, the methodology can be applied to the examina- 
tion of nonlinear waves [26]. Since this entire topic is very broad and still very active, a com- 
plete survey is not within the scope of this report. Rather an outline of some of the features 
of the various methods will be undertaken, and a more detailed presentation of the methods 
most familiar to the author will be presented. Only linear waves will be considered here. 

Two basic approaches to the study of waves in anisotropic and dispersive media have 
been developed. The first approach has been called the study of the kinematics of waves [24, 26] 
and also the description of waves in terms of rays [4]. The approach is very powerful in its most 
general form when local properties of the waves are sought. When boundaries must be 
considered, it is not clear that the method can be applied as conveniently. It is this theory, 
when pushed to greater generality, which provides methodology for examination of nonlinear 
waves. An outline of the simpler aspects of this approach will be presented in this subsection. 

The second approach is a Fourier analysis approach and has been called the method of 
modes by Turner [4]. It has been developed to a highly useful form by Lighthill [25] and by 
others. In this subsection the method of Fourier analysis, including extraction of relevent 
information by asymptotic evaluation, will be emphasized because this general methodology 
has been found so useful in applications. Following this subsection, some applications will be 
presented. 

a. Kinematics of Waves 

As pointed out by Lighthill [18] and discussed in great detail by Whitham [26], the method 
of the kinematics of waves "makes use of the duality between waves and particles, that is 
fundamental to quantum mechanics, and follows a wave packet through the system as if it 
were a particle." The equations for the waves are developed in Hamiltonian form, a Hamilton- 
Jacobi equation is found for the phase of the wave, a Lagrangian is found for the wave motion 
and a variational approach to the determination of the equation is developed. One advantage 
to this approach is that the P.D.E. for the waves can be cast in the form of ordinary differential 
equations (O.D.E.) which are much easier to solve computationally. A second advantage is 
that this approach provides a systematic methodology for attacking nonlinear problems. 

In the one-dimensional case, a statement of the conservation of waves in differential form 
can be written as 

S+!-°- < 3 °> 

This equation states that the time rate of change of the number of waves per unit length (the 
wave number) plus the gradient of the flux of waves (the frequency) out of an element equals 
zero. If, in addition, a dispersion relation is known for the wave system, which in general depends 
upon position and time as well as wave number, then 

a=W(k, x, t) (31) 

and the continuity equation can be written as 

dk./dW\ dk_ -dW , Q9 . 

dt + \dk)di~~dx~' [6Z) 



In this form the equation is hyperbolic with a characteristic defined by 

dx dW 
and characteristic form 



dt-» (33) 



dk_ dW m) 

dt w (34) 

These two ordinary differential equations, obtained by writing the kinematic conservation 
equation and the disperson relation in characteristic form, are simply the equations of motion 
for a particle written in Hamiltonian form. The following identifications can be made : x is the 
position of the wave in one case and of the particle in the other; just as in quantum mechanics, 
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the wave number k is identified with momentum p and the dispersion function W is identified 
with the Hamiltonian H. Furthermore, the frequency co of the wave packet serves the same role 
as the energy E in Hamiltonian mechanics, so that 

dt~ dt ~ dk dt + dx dt + dt 

do> = _dW dW dW dW.dW 

dt dk dx + dx dk + dt {6b) 

do) 3co 

di^m' 

From these equations it can be seen that the wave packet travels with the group velocity 

_dw 

Cg dk' 
When the medium is homogeneous, the dispersion relation is independent of position and 

§=<>• 

Hence the wave number (momentum) is constant along a characteristic (trajectory). When the 
dispersion relation is time independent, then 

-"-=0 
dt U 

and the frequency (energy) is constant along a characteristic. 

One more identification between waves and Hamiltonian mechanics will be made. The 
frequency and the wave number can be defined as partial derivatives of the phase function 
6 as follows 

w= -§r k= ¥x m 

Here 9 can then be identified with the action in mechanics. With these definitions, the dispersion 
relation becomes the Hamilton-Jacobi equation when w and k are eliminated : 



w=W(k,x, t) 
becomes 

dd 
dt 



+^(|>m)=0. (37) 



All of these equations can be generalized in a straightforward fashion to the three-dimen- 
sional case, and these equations will simply be listed. 

w=W(k, x, t) (38) 



or, in component form 



dx 1= dW 
dt dki 

dk 1= _dW 
dt dx l 
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dx 2 
It' 


dW 

dk 2 




dk 2 _ 
dt~ 


_dW 
dx 2 




dx 3 


dW 




dt 


~dh 




dk z _ 
dt 


_dW 

dk 3 




dd 
u= di' 


k=v x e 


de 
dt 


4-W(v*0,M)=O 



(39) 



(40) 

and 

an 

(41) 

In this form (O.D.E.), these equations are ideally suited for ray-tracing techniques using com- 
putational methods. 

Finally, as noted before, this formulation of waves can be extended by considering the 
Lagrangian for the system and using variational methods with the Lagrangian. This method- 
ology is very well suited for examination of nonlinear waves. Since it is an area of active research 
and application and since the scope of this paper is limited, the interested reader is referred 
to reference [26] for the extensions. 

b. Fourier Analysis and Asymptotics 

In this subsection, the outline of a Fourier analysis of a class of linear P.D.E.'s describing 
waves in anisotropic, dispersive media will be given. The analysis was presented in this form by 
Ligh thill [25], although similar analyses, of somewhat less generality, have been presented 
elsewhere. The analysis, as presented, is appropriate only to infinite homogeneous media, but 
extensions to inhomogeneous media or media with boundaries will be mentioned at the end. 
Applications of this analysis will be discussed in the following subsections. 

For a general infinite, homogeneous medium in an undisturbed state, a small disturbance 
is governed by a linear P.D.E. with constant coefficients, which may be written 

K4-«;r-w- ; l)*=° w 

where P is a polynomial and is a dependent variable specifying the wave. A plane-wave 
solution 

0=0o exp {i(— ut+k-x)} (43) 

can exist if the dispersion relation 

P(^k u k 2 ,k 3 )=0 (44) 

is satisfied. 

When the disturbance is produced by a localized forcing region, then the right side of 
eq (42) must be replaced by a "forcing term." For this inhomogeneous form of the equation, 
it is necessary to superpose the various plane-wave modes for the homogeneous equation to 
obtain the solution. A very general form for such a forcing function is 

exp(-W)/(x-U0 (45) 

where co is the frequency with which the forcing term oscillates and U is the velocity with which 
it travels. This forcing term represents one which translates and oscillates simultaneously 
exciting waves in the otherwise undisturbed medium. A slight modification to the functional 
form above allows one to consider transient phenomena as well as steadily time varying dis- 
turbances [25]. Here only forcing functions of the form above will be considered. 
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By a localized forcing function it is meant that /(x) =f(x, y, z) vanishes outside some region 
around the origin. Then it can be decomposed into a Fourier integral 

/( X )= r dh r dh r ^ 3J w k - x (46> 

J — 00 J — co J — oo 

where F(k)=F(k lf k 2 , h) is a regular function of k. The equation 
has the formal solution, using the Fourier-integral procedure, 

-;:.*/:*/:* m> ^wr vmi 

This solution to the inhomogeneous equation is not unique: any solution to the homo- 
geneous equation can also be superposed to produce another solution. However, when the 
requirement is imposed that only outgoing waves appear at infinity, the so-called radiation 
condition, then the solution becomes unique. Lighthill [17, 25] has given a prescription for 
applying the radiation condition and making the solution unique. The denominator of the 
integrand in eq (48) is the dispersion relation for waves excited by the forcing function. Points 
at which this function vanishes arc poles of the transform kernel. Using residue theory, one of 
the Fourier transforms can be inverted exactly, and the radiation condition can be applied 
during this inversion [17, 25]. 

The solution, after one inversion has been performed, is still very formidable. Under 
special circumstances the other two integrals can also be inverted in terms of known functions; 
however, such inversions are the exception rather than the rule. As a result, asymptotic evalua- 
tion of the solution in the far field is very important for examination of the wave system. 
Again, Lighthill [17, 25] has given a prescription for such asymptotic evaluation using a multi- 
dimensional (two-dimensional in this case) stationary-phase analysis. 

Since F(k) is a regular function of each k u F(k) has no singularities and can be taken as a 
smooth, well-behaved function; therefore it remains passive in the asymptotic analysis. The 
key to the asymptotic analysis is the dispersion relation 

P(co +U.k, k)=0. (49) 

For fixed co and U, the dispersion relation describes a surface in wave number space. At regular 
points on this surface, points at which the curvature of the surface does not vanish, a two- 
dimensional steepest-descents procedure can be applied. This procedure determines both the 
directions in real space into which waves spread and the rate of decay of such waves with dis- 
tance from the forcing function. 

When the curvature of the dispersion-relation surface vanishes, a more involved steepest- 
descent procedure must be used. In the directions in real space for which this occurs, the ampli- 
tude of the waves fall off more slowly with distance from the forcing function than in other 
directions. For details of the asymptotic procedure, references [17] and [25] should be consulted. 

When inhomogeneous media (media whose properties depend upon position) or finite 
media (media with boundaries) are considered, the problems become more difficult and often 
recourse must be made to numerical methods. Under these conditions, the inhomogeneity or 
the boundaries can change the spectrum of the eigenvalues from continuous to discrete. In the 
next subsection on stratified-fiuids, where boundaries and inhomogeneities are important, this 
will be discussed in more detail. 

3.3. Waves in Stratified Fluids 

In this subsection a surve}^ of internal waves in stratified media will be presented. The 
emphasis will be upon waves in a quiescent, stratified fluid using a modal analysis, as discussed 
in the previous subsection. First, the equations will be derived. Then a discussion of results 
will be presented. 

a. Derivation of the Linearized Equation for the Vertical Displacement 

The equations for mass and momentum conservation are eqs (2.11) and (2.12). [Note that 
equations will be numbered consecutively in each section of this report. When equations from 
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a preceding section are referenced, the equation number will be prefixed by the section number, 
as the equations referenced above.] Equation (2.11) will be written as 

^+pV'U=pQ. (50) 

Here Q is a source term; this term has been included to discuss the modeling of internal waves 
(IW) excited by body motions or by "wake-collapse." Equation (2.12) will be written in a 
nonrotating reference frame 

jjf+lvp-9i>=0 (51) 

where g = — gi z and i 2 is a unit vector in the z direction, which is taken to be positive upward. 
The equation of energy conservation, eq (2.16), will be rewritten with the energy addition 
term set to zero. This isentropic relation then becomes 

!+u.V^(^+u.Vp) (52) 

For an incompressible fluid the speed of sound c becomes infinite, or the density following a 
fluid element does not change. Hence 

and from the continuity equation 

V-u-Q. (54) 

In equilibrium the fluid is quiescent and stratified in the vertical direction, so that 

Uo=(w , v , w )=0 

Qo=o 

(55) 

Po=Po(z) 



Po(z)— p (z=0) =—g J p Q {z')dz f . 



The disturbance quantities are denoted by primes, 

Q=Q' 

(56) 
p=Po(s) + p / 

P=Po(z)+p'- 

These quantities are substituted into eqs (51), (53), and (54); and the equations are linearized 
with respect to disturbance quantities. The prime notation on disturbance quantities is dropped, 
and differentiation is denoted by subscripts. Then eqs (54), (53), and (51) become 

u x +v y +w z =Q 

Pt+wpo,z=0 (57) 

1 1 1 gp 

U t = Px, V t = p V) W t = Pz~— ' 

Po Po Po Po 
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These equations can be reduced to a single equation for the vertical displacement f of a 
fluid element from its equilibrium position : 



In this equation 
and 



(d t 2 +m (r„+jw)+r M „=G,,+^ (&-?.„). (58) 

Po 

N 2 =-g( Po>z / Po ) 



is the square of the Brunt-Vaisala frequency, which, in general, is a function of the depth z. 
The last term on the right is taken to be zero by the Boussinesq approximation; therefore the 
final equation for the vertical displacement is 

(6y+Af 2 ) (r«+r w )+r„„=&,. (59) 

Equation (59) agrees with the linear portion of eq (5.24) in reference [5]. Some authors have used 
other dependent variables or pairs of equations to formulate problems in stratified fluids. 

To complete the formulation of problems in stratified fluids, boundary conditions must be 
applied. If the fluid is considered to be infinite and the source of waves localized, then a radiation 
condition (only outward propagating waves) must be applied at infinity. If the fluid is bounded 
by a solid surface, then the velocity normal to the surface must vanish. If there is a free surface, 
then the pressure remains continuous across the free surface. 

Finally, initial conditions must also be specified. Since eq (59) is second order in time, 
both the initial displacement field and the initial vertical velocity field can be specified. The 
initial conditions can be applied formally if Laplace transform techniques are utilized. If 
steady state waves only are of interest, the initial conditions and transient effects may be 
ignored. 

b. A Survey of Results on Internal Waves in a Quiescent Stratified Fluid 

In an incompressible fluid the stratification can vary in a variety of ways. The simplest 
case to anatyze is that in which the Brunt-Vaisala frequency is constant and the fluid is infinite. 
(A constant Brunt-Vaisala frequency implies that the fluid is exponentially stratified with 
density 



Po(2)=Poo expf — -z\ 



where p 00 is the density at 2 = and z is measured positive upward.) This case also illustrates 
many of the interesting phenomena which occur in stratified fluids. In this subsection a survey 
will be given of some of the books and papers which describe internal waves; this survey will 
include a summary of methods and results, but will not include mathematical detail. 

As shown above, eq (59) governs the internal-wave motion in a Boussinesq fluid. For an 
infinite fluid the plane-wave solution will satisfy this equation: 

f=Aexp {i(-ut+ax+(3ij+yz)}. (60) 

(The notation has been changed from that used in subsection 3b; a, /3, and y are used as wave 
number components rather than k h k 2 , and & 3 .) This plane wave solution satisfies eq (59) for 
nonzero A if the dispersion relation is satisfied: 

(a 2 +P 2 )(N 2 -a> 2 )-yW=0. (61) 

From this relation it can be seen that plane waves, with a, 13, and y all real, only exist provided 
that u<N. If co>iV, one of the wave numbers must be imaginary for the dispersion relation 
to be satisfied, and internal waves cannot propagate. Therefore, iV is a high-frequency cutoff 
for internal waves. 

If the wave numbers are written in spherical polar coordinates 

a=k sin 6 cos <£>, f3=k sin sin <£, y==k cos (62) 
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where k is the magnitude of the wave number and 6 and <i> are the angles defining its direction, 
then the dispersion relation becomes 

tan 2 6= 1 tti 5 



or (63) 

sin 2 0=co 2 /iV 2 

Therefore the frequency of the wave does not depend upon the magnitude of the wave number, 
but only upon its direction. The anisotropic nature of the medium can be clearly seen from 
this form of the dispersion relation: for a plane wave of fixed frequency the direction of the 
wave number with respect to the vertical axis in wave number space is fixed. 
The phase velocity, defined in eq (24), is 



\a 

']k\ 2tx ~(« 2 +P 2 +y 2 ) 



c= ^k=7ICTT-^ \i\ ( 64 > 



TJ 



and the group velocity, defined in eq (29), is 



1 (N 2 -co 2 ) 2 



Y «' 



(« 2 +£ 2 ) 



(65) 



The dot product of c and c g shows the rather unusual property that the phase velocity and the 
group velocity of these internal waves are orthogonal. Therefore, energy is propagated at right 
angles to the velocity of an individual wave. This property is a result of the fact that the medium 
is both anisotropic and dispersive. 

Gortler [27] examined the steady-state wave configuration in a homogeneous medium of 
constant Brunt-Vaisala frequency. He discussed the properties of the wave system in terms 
of characteristics of the linearized equations of motion noting that N is the high-frequency 
cutoff for internal waves. Using a Schlieren system, he determined the wave-configuration 
experimentally. 

In an independent study, Mowbray and Rarity [28] examined the two-dimensional internal 
waves excited by a cylinder oscillating in a constant- N stratified medium. In their experiments 
they duplicated the results of Gortler using a schlieren system and explained these results in 
terms of the energy propagated by the group velocity in the anisotropic, dispersive medium. 
They found that, when internal waves are excited, significant amplitudes occur only in regions 
emanating along preferred directions from the oscillating cylinder. These directions are deter- 
mined by the angle 6, measured from horizontal, where 

smd=±u/N. (66) 

For a specified frequency co, 6 is the angle of the group velocity. The variation of the angle d 
with variation of frequency co was measured, and N was verified as the high-frequency cutoff. 

From the form of the plane-wave solution (60) and the dispersion relation (61), 
it is apparent that the eigenvalues of this system form a continuous spectrum. From the prin- 
ciple of superposition, an integral over the eigenfunctions (60) form the general solution to 
problems in an infinite, constant-iV medium, as discussed previously. When boundaries are 
present or the Brunt-Vaisala frequency varies with depth, the eigenvalues can change from a 
continuous to a discrete spectrum. 

To illustrate this behavior, a layer of fluid of finite depth D will be considered. In this 
case, the general solution can be written as a function of z times a horizontal plane wave: 

f=Z(s) exp {H-ut+ax+Py)} (67) 

and the homogeneous version of eq (59) becomes 

oo 2 Z" + [N 2 (z)-u 2 ] (a 2 +/3 2 )Z=0 (68) 
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where the prime denotes differentiation with respect to z. This equation corresponds to eq 

1 y 

(5.2.8) in reference [5]. The boundary conditions for eq (68) are that the vertical velocity -j- 

vanish at z=—D and that the pressure remain constant (at atmospheric pressure) at the free 
surface 2=0. Phillips [5] presents the proper boundary condition to first order at the free surface; 
but for internal waves this surface can be treated approximately as a rigid boundary [29]. 
Therefore, at z = 0, the vertical velocity vanishes, and, in terms of the vertical displacement, 
£=Z=0 at z=0 and at z=—D. 

For general variations of A (2), solutions can be found only for discrete values of the fre- 
quency o) n . Corresponding to each eigenvalue co n is an eigenfunction \p n (z), and any solution to 
the inhomogeneous form of eq (68) can be expanded in terms of these discrete eigenf unctions. 
If N is constant in the layer for example, 

. nwz 
\p n =sm -jj 

and 

11/2 

/ TVIU \ o , „ 

N' 



L « 2 +/3 2 



2 + /5 2 



An inhomogeneous form of eq (68) arises when disturbances are produced by nonzero initial 
conditions on the displacement or its time derivative, or by some internal forcing function. 

Several authors have examined waves in inhomogeneous, stratified media. Eckart [6] has 
discussed the general methods for analysis presented in earlier subsections: he has discussed 
both ray-tracing techniques and modal techniques for analyzing waves. Also, he has discussed 
the effects upon the wave systems of compressibility and of spherical geometry. Furthermore, 
he has considered inertial waves in rotating fluid and some of the aspects of the infinitesimal 
waves in rotating and stratified fluids. Reference 6 is the most complete examination of waves 
in stratified and/or rotating fluids with application to geophysical fluid dynamics. No effects 
of turbulence or viscosity are considered. 

In a chapter on waves of small amplitude, Yih [7] has examined stratified, compressible 
fluids, concentrating on fluids of finite depth. He has examined several mathematical questions 
related to problems that arise in studies of these fluids. 

Phillips [5] has discussed the internal waves in a finite layer of incompressible fluid when 
the variation of the Brunt-Vaisala frequency models a thermocline. He derives eq (68) and notes 
some of the properties of this equation for a peaked distribution of N(z). Equation (68) is in a 
form such that Sturm-Liouville theory can be applied directly; the modes of this equation which 
can be excited depend upon the value of co with respect to N max , the maximum value of the Brunt- 
Vaisala frequency. When co>A^ max , no internal waves are excited, and this is a statement that 
Amax is the high-frequency cutoff. When co<^V max , discrete values of co exist as eigenvalues for 
the problem. In the lowest of these modes for internal waves, Z(z) has its maximum near the 
maximum of N(z), that is, near the thermocline. In this mode the whole thermocline heaves up 
and down in phase. The second eigenfunction or mode undergoes one change of sign in the ther- 
mocline region. This mode represents an internal wave which thickens and thins the thermocline 
region. Each higher mode has one more change of sign with oscillatory behavior in between in 
the thermocline region. 

Several other aspects of internal-wave phenomena will be discussed briefly. Excitation of 
internal waves, waves by a body traveling through a stratified fluid, has been an area of rather 
active applied research during the past few years. It has been found that there are three ways 
in which a moving body can excite internal waves in an incompressible fluid. First, because of 
the body motion, a fluid element directly in the path of the body will be displaced as the body 
passes, and this displacement will produce internal waves. This effect has been studied both 
theoretically and experimentally by many people. 

A key study was performed by Miles [29] who calculated the magnitude of the internal 
waves generated by this effect in both a constant- JV model and a thin-thermocline model of the 
fluid. Important features of this analysis are that the internal waves are related to motion and 
the size of the body using a small, slender-body approximation and that a stationary-phase 
analysis is used to calculate the far-field properties. Carrier and Chen [30] analyzed this effect 
in the far field for a three-layer model for the thermocline structure and for a finite size body. 
Experimentally, the internal-wave phase configuration was observed first by Stevenson [31]. 
Additional work in this area has been concerned with refining theoretical models of the thermo- 
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cline structure, examining the interaction of internal waves with a fluid surface and with 
making measurements of internal-wave amplitudes. 

A second mechanism by which a moving body may produce internal waves is by " wake- 
collapse" [32]. Any body moving through a fluid produces a turbulent wake due to frictional 
drag and, if it propels itself, due to the propellor action. The turbulent wake grows with dis- 
tance behind the body, mixing the ambient fluid into the wake and tending to homogenize the 
fluid. At a certain distance downstream from the body, growth of the wake ceases, and, due 
to the buoyancy effects, the wake begins to collapse. During collapse, the fluid in the wake 
moves to the level at which its density is equal to the ambient density. The collapse is a 
reasonably efficient generator of internal waves. 

Studies of wake collapse and internal-wave generation have proceeded along two lines. 
One line of study is the examination of the turbulent-wake properties and scaling of these 
properties with body size, speed etc. [33, 34, 4]. 

The other line of study has emphasized the internal waves produced by such collapse 
[35, 29, 30]. Most of the work along both lines has appeared in report form only; the references 
given here are indicative of the type of studies being carried out. 

A third mechanism by which internal waves may be excited by a body in a stratified fluid 
is through oscillations of the body at a frequency below the Brunt-Vaisala frequency [28]. 
When the body remains at a fixed location and bobs up and down or "breathes fluid' ' in and 
out, these waves have been analyzed [36, 37] and measured experimentally [38]. When the 
body is both translating and oscillating at a prescribed frequency, either greater or less than 
the Brunt-Vaisala frequency, the internal-wave system has been analyzed [39] and the phase 
configuration of the waves have been observed experimentally [40]. 

Limitations of space do not allow any discussion of other currently interesting areas of 
research concerning small-amplitude internal waves. The books of Turner [4] and Phillips [5] 
discuss some of these, for example, internal waves in moving stratified media and weak nonlinear 
interactions between waves. 

3.4. Waves in Rotating Fluids 

In this subsection, a survey of inertial waves will be presented. It will be rather brief 
because many of the methods for analysis of the waves were discussed before and because many 
of the results are the same as those in stratified fluids. First, the equations governing small- 
amplitude waves in rotating fluid will be derived, and features of intertial waves in an uncon- 
fined, uniformly rotating fluid will be presented. Then a short discussion of Rossby waves, 
which arise in more general rotating flows, will be given. 

a. Derivation of Equations 

From the equations for mass and momentum conservation, eqs (2.11) and (2.12), the basic 
equation for the vorticity of a fluid element can be derived [9, 10]. Formally, the vorticity is 
defined as the curl of the velocity vector, and physically it is twice the angular velocity of 
rotation of a fluid element. Equation (2.12) can be rewritten as 

^+\ Vu 2 -uX (a,+2!2)= — v^+V [-03+^ (J2Xr) • (0Xr)l (69) 

When the curl of this equation is taken, an equation for the vorticity w =vXu results 

a 

Using eq (2.11), assuming 



a "-VX[uX(cu+2I2)] = -v(-)xv^. (70) 



^=0 
dt u 

and manipulating somewhat, eq (70) can be rewritten as 

Equation (71) is an equation for the absolute vorticity [9, 10], w + 2Q, which is the vorticity 
in an inertial (nonrotating) reference frame. It states that the time rate of change of the absolute 
vorticity divided by the density following a fluid element is influenced by two processes. The 
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first is simply the change in vorticity due to deformation of the vortex lines (stretching and 
tilting) as the vortex lines are convected about with the fluid (see ref. [9] and references therein 
concerning vorticity). The second term is a production term for vorticity resulting when the 
gradient of the density and temperature are not parallel. When these gradients are parallel, 
this term is zero, and pressure is a function of density. (When viscous effects are also important, 
another term is present on the right-hand side of eq (71), and this term is diffusive in nature. 
Several standard references on fluid dynamics have discussions of the viscous diffusion of 
vorticity) . 

The vorticity is a very useful concept, and eq (71) is an important equation from which 
many useful relations can be derived. For example, the Taylor-Proudman theorem, (Eq 2.13a), 
follows directly from this equation when the flow is incompressible (or barotropic so that 
pressure is a function only of density), inviscid, steady and slow (so that products of vorticity 
co and velocity u can be neglected). 

Equation (71) can also be used to derive other governing relations, such as the fundamental 
equation for small-amplitude, inviscid inertial waves in an unbounded, incompressible, rotating 
fluid. In an incompressible fluid V-w = 0, and, therefore, by the continuity equation (2.11), 

Dt U 

Then, formally the density can be removed from within the substantial derivative in eq (71). 
Since 12 is constant, its substantial derivative is zero, and because nonlinear terms are being 
neglected, the substantial derivative of the vorticity reduces to the partial of this quantity 
with respect to time. On the right-hand side, the last term is zero because the fluid is assumed 
to be barotropic, and in the first term the nonlinear portion w-vw is neglected. Therefore, 
eq (71) becomes 

^=2fi.vu. (72) 

dt 

Taking the curl of eq (72) and eliminating the vorticity between the new equation and eq (72) 
yields a single for the velocity 

[^ 2 V 2 +(2fi-v) 2 ]u=0- (73) 

This equation is the fundamental one for inertial waves in an inviscid, incompressible, 
uniformly rotating fluid [10]. As in the preceding subsection, an excitation source could have 
been included in this derivation, resulting in an inhomogeneous form for the equation. However, 
discussion will be restricted to the homogeneous form. The formulation is complete when initial 
and boundary conditions are specified. For an infinite fluid and steady-state considerations, 
only a radiation condition needs to be applied at infinity. 

Another important and interesting equation can be derived from eq (71), an equation 
which governs the vertical component of vorticity in a rotating layer of incompressible fluid. 
In this approximation two key ideas are utilized. First since the fluid is rotating, the absolute 
vorticity obeys eq (71) (with the second term on the right-hand side zero) and the motion is 
essentially two dimensional, as the Taylor-Proudman theorem states when the motion is slow 

A A 

and steady. If fl=S2k where k is a unit vector in the positive z direction, then the motion is 
such that the fluid variables approximately only depend on x and y. Then columns of fluid 
move about as units such that the horizontal components of velocit}^ are uniform with height. 
Second, since the fluid is incompressible, changes in the horizontal divergence of the fluid 
are accompanied by changes in the fluid height. Then the basic unit to consider to describe 
the fluid motion is the fluid column. If H(x, y, t) is this height of the fluid, then the incompres- 
sibility condition v-u = becomes 

/du ,dv\_dw_ 1 DH (7 , 

\dx + dy/dz~HlM { } 

where 

Dt 
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is now the two-dimensional substantial derivative 

Dt~dr U dx^ V dy 

A 

With f=w-fc, the vertical component of eq (71) becomes 

^ (tt-20) = (H-20) ~ (75) 

where, again, the substantial derivative is two dimensional since, as discussed above, u and v 
have been assumed to remain approximately constant with height. 
Combination of eqs (74) and (75) yields 

The quantity (f +20)/£Tis called the potential vorticity [9, 10, 73, 74], and the relation (76) 
is known as the conservation of potential vorticity. It states that the z component of the abso- 
lute vorticity divided by the fluid layer thickness is constant following a fluid element (column). 
Its applicability to the fluid layers on a spherical surface for representing the oceans or atmos- 
phere, is a topic of considerable interest [9, 10]. 

In the next subsection a brief discussion of inertial waves, starting from eq (73), will be 
presented. Then in the following subsection, eq (76) will be examined to show how another 
type of wave, Rossby waves, arise. 

b. A Survey of Results on Inertial Waves in Unconfined Uniformly Rotating Fluids 

For an infinite fluid the plane wave solution 

U=V exp [i(— tot+ax+Py+yz)] (77) 

will satisfy eq (73) provided that the dispersion relation 

» 2 Jfc 2 -4QV=0 (78) 

where 

Jc=(a 2 +f3 2 +y 2 y /2 

is the magnitude of the wave number. Equation (78) shows that solutions exist for real a, fi, 
and y only for co<212, exactly as in the stratified-fluid case with N replaced by 212. If spherical 
polar coordinates are used in wave number space, then y=k cos 6 and 

cos 2 B=~- (79) 

Comparison of eqs (79) and (63) shows that the dispersion relations are the same with 212 
replaced by N and cos 6 replaced by sin 6. Therefore, the observations made about internal waves 
apply also to inertial waves. For example, the phase velocity is 

fa) 

(80) 



(a 2 +/3 2 + 7 2 ) |;j 
and the group velocity is 

2 412 V \ ~ a \ /oin 

co (a +P+y) [( a 2 + p2)/ y \ 

Therefore, as in the case of internal waves, the group velocity is orthogonal to the phase velocity 
for inertial waves. 

The first discussions of the steady wave system in a rotating fluid were carried out by 
Gortler [41], who noted that the characteristic directions for the waves are given by the relation 

cos0=±co/212 
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where is the angle in physical space measured from horizontal. Subsequent analysis of the 
steady-state inertial waves excited by an oscillating disk in a rotating fluid were made by Oser 
[42], who then verified the analysis experimentally [43]. In his analysis Oser was able to express 
the solution for the small-amplitude waves in terms of known functions and to calculate the 
flow field near the disk as well as in the far field. 

For a fixed rotational frequency ft, the oscillation frequency co determines the direction of 
energy propagation for the inertial waves, eq (82). When the oscillation frequency is small, 
this direction becomes nearly vertical, whereas when co— >212, the direction is nearly horizontal. 
This behavior is to be contrasted with that in a stratified fluid. There the characteristic direc- 
tions become nearly horizontal for small co and nearly vertical as co— >N. 

No additional studies on waves in rotating fluids will be discussed here due to space limi- 
tations. Some additional considerations of inertial waves in unbounded media are given by 
Greenspan [9], who, discusses, among other questions, the reflection of plane inertial waves 
by an interior surface, the slow transient motion of a disk along the axis of rotation (and how 
the phenomenon known as a Taylor column develops) and waves produced by a travelling 
forcing function (along the lines presented earlier for a stratified medium). Eckart [6] also 
presents analyses and results on waves in rotating fluids. Furthermore, he discusses small- 
amplitude waves in fluids when rotation and stratification are both present. Some more recent 
studies of waves in rotating and stratified fluids have been performed by Hendershott [36], 
who examined waves excited by an oscillating body, and by Rao [44], who examined waves 
excited by travelling forcing effects. 

c. Rossby Waves 

There exists in large scale geophysical flows another type of wave, the Rossby wave, which 
owes its existence to a mechanism other than that described in the last subsection. Here an 
attempt will be made to explain this mechanism, and some of the features of these waves will 
be discussed. More detailed discussions of Rossby waves, with consideration of laboratory 
experiments and of oceanographic and atmospheric waves, are given in references [9, 10]. 

The conservation of potential vorticity, eq (76), is a statement that the component of the 
absolute vorticity in the direction of the rotation vector divided by the depth of a fluid column 
is constant to an observer following the column. On the earth, the direction of the rotation 
vector is from the south to north poles so that the conservation of potential vorticity is strictly 
valid only in this direction. Often the conservation relation is applied for the component of 
absolute vorticity in the radial direction at that latitude. When this is done, 212 in eq (76) is 
replaced by/=2£2 sin 0, where 6 is latitude, and/ may be either considered to be locally constant 
or to vary with latitude. The conditions under which the conservation of vorticity applies in 
the radial direction have been discussed in reference [9] for example. 

Rossby waves can arise when either/ or H varies with lateral position. For example, for 
relatively large scales on the earth,/ varies approximately linearly with latitude; in the northern 
hemisphere, / will increase in the northerly direction (as 6 increases) : 

f&fo+PV (83) 

where y is taken to increase to the north and x to increase to the east. Substitution into eq (76) 
with H constant, and linearizing gives 

A similar situation occurs when the imposed angular velocity 212 is constant and the depth of 
the layer varies linearly with external position [9]. Again, eq (84) is approximately the governing 
equation, where /3 is now a measure of this linear variation. 

Equation (84) is simply a linearized statement of the conservation of potential vorticity. 
The variation with lateral position of the imposed angular velocity, eq (83), provides a restoring 
mechanism, and it is this restoring mechanism that gives rise to Rossby waves. For any dis- 
turbance, the rate of change of position in the northerly direction, v } produces a corresponding 
reduction of the time rate of change of its vorticity because of the variation in the imposed 
angular velocity/ 

Equation (84) can be rewritten in terms of a single dependent variable, the stream function 
\f/\ where 

um *k, vm *k. (85) 

dy dx 

215-653 O - 77 - 4 
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Then 



and eq (84) becomes 



dv du 2 . 



|*H*g-0 (86) 

a form examined in some detail by Lighthill, using Fourier analysis and asymptotic techniques. 

4. The Analogy Between Stratified and Rotating Fluids 

As discussed in the introduction, there is an analogy between stratified and rotating flows. 
Stratification and rotation provide constraints on the fluid flow, the stratification inhibiting 
the flow in the vertical direction and the rotation inhibiting flow in horizontal planes. In section 
3 on Waves, this analogy was given mathematical form. There it was shown that internal waves 
in a stratified fluid of constant BV frequency behave like inertial waves in a rotating fluid. 

In this section the mathematical analogy demonstrated in the previous section is extended. 
The review article of Veronis [16] on this analogy is used as the basis for this section. Since that 
article is both complete and clear, this section will be relatively brief and will concentrate only 
on the derivation of the equations from which the mathematical analogy is identified. 

The motivation for examination of this analogy is provided by the desire to have a deeper 
understanding of both physical phenomena. Understanding phenomena in one system can 
clarify the analogous phenomena in the other system. In addition, laboratory studies in one 
system may be much more easily carried out than in the other system, allowing problems to be 
studied in the simpler system. 

The systems considered will be assumed to have the axes of rotation and stratification 
coincident in the vertical direction. Dissipative effects such as viscosity and thermal conduc- 
tivity will be included in the analyses. Also, boundaries are considered, but only those which 
are parallel or perpendicular to the rotational and gravitational direction. The analogy will also 
be seen to hold for systems which are both rotating and stratified. 

In section 3, it was noted that waves in stratified and rotating fluids are related when the 
identification between the Brunt-Vaisala frequency N and twice the rotational frequency 212 
is made. Closer examination of the perturbation equations shows the further identifications 
which must be made in analogous flows. Equations (3.57) are the perturbation equations for 
waves in a stratified fluid. If it is assumed that there is no dependence of the waves upon spatial 
variables and that no sources are present, Q = 0, then these equations become 

u t =v t =0 (1) 

PoW t =—gp. 

Here, the continuity equation is identically satisfied, p (z) is the equilibrium density distribu- 
tion, p is the perturbation density and u, v, and w are the perturbation velocity components. 
From these equations, the velocities u and v in the horizontal direction are seen to be constant 
and taken to be zero. Only the remaining two equations for the density perturbation and the 
vertical velocity perturbation are relevant: 



Po Tt = ~9P 
dp 

The analogous perturbation equations can be obtained for a rotating fluid. If it is again 
assumed that there is no dependence upon spatial coordinates, then in component form these 
equations are 

at 
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%=- 2 ® u (3) 

dt u - 

The last equation implies that the vertical velocity is constant, and it is taken to be zero. 
Only the first two equations are relevant. 

Equations (2) and the first two of eqs (3) have exactly the same form; in fact, they can 
be made exactly the same if, in eqs (2), u is substituted for w, v is substituted for (—(Jp/pqN) 
and 2Q is substituted for N. From eqs (2) and (3) and other simple examples of this nature, 
the following identifications, except for multiplicative factors, can be made [16]: 

Rotating s}^stem (fl) v w u z x 
Stratified system (S) —p u w x z. 

The more general equations expressing the analogy for fluids which are both stratified and 
rotating will be derived following reference 16. It is assumed that the stratified fluid can be 
described in the Boussinesq approximation, and other restrictions stated above also apply. 
The equation for the conservation of mass (2.11) becomes 

V-u=0. (4) 

The equation for the conservation of momentum (2.12) becomes 

^+u.Vu+2i2Xu=--vP+^TL+^V 2 u. (5) 

Ot po 

Here P is the reduced pressure (including the centrifugal-force term and the static part of the 
gravitational-force term [9, 16]). The buoyancy term is taken in a form appropriate for a 
Boussinesq fluid: the density is taken to be a linear function of the temperature variation T 
from a reference value (see eq (2.24) and the discussion of that equation). The force f in eq (2.12) 
is taken to be the viscous force *>V 2 u. 

The equation for the conservation of energy is eq (2.19). If the temperature is written 
as T-\-T, where T(z) is the temperature distribution in the absence of motions, then eq (2.19) 
becomes 

^+u^T+w^=kv 2 T. (6) 

dt dz 

Here it has been assumed that T(z) is a linear function of depth and that T is the perturbation 
temperature caused by the motion. Equation (5) already implies this separation of the tem- 
perature. The constant temperature gradient will be defined to be equal to 

rfT =4 AT 
dz L 

AT is the magnitude of the temperature difference over the vertical length scale L, the factor 4 
being included for convenience. Then eq (6) becomes 

^+u-vT+±^w=kV 2 T. (7) 

Equations (4), (5), and (7) will be nondimensionalized in two ways, depending upon the 
physical system to be modeled. For the first system, denoted Sq, stratification is dominant 
with rotational effects modifying the motion. For the second system, denoted fl s , rotation 
dominates and stratification modifies the basic motion. 

A general nondimensionalization for both systems can be made using the following scales 
and nondimensional variables: 

u=y U / , t = T t', r=Lr', P=PP' } T=26T'. 
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Here, the primed variables are dimensionless, and V, r, L, P and 6 are the dimensional scaling 
factors. Equations (4), (5), and (7) become 

V'-u'=0 (8) 

7^ + T u ' -v ' u ' +2 " FizXu ' =_ iz V'J"+20&07"i,+ J W (9) 

^ + ^ v 'T'+^ W '=f ^ T , (10) 

4.1. The Stratification-Dominated System, Sq 

When stratification dominates, the various dimensional scaling factors are related by 
equating the magnitude of the time acceleration term, the pressure term and the buoyancy 
term in eq (9) Therefore, 

V P 

T-JXT** (11) 

Also, the magnitude of the time derivative of the temperature and the convective temperature 
changes due to the equilibrium temperature distribution are equated in eq (10). Then 



26 2ATV 
r~ L 

Equations (11) and (12) imply the following relations between scaling factors: 



(12) 



(13) 



d=^AT, P= Po Lgae. 

Equations (8)-(10) become, after dropping the prime notation, 

V-u=0 

^+eU-VU+ji 2 XVL=-vP+2Ti 2 +a 1/2 Rv 2 U (14) 



where 



^+eU'VT+2w=a~ 1/2 Rv 2 T 
ot 

€= X = aT F== sTt (15) 



a =vjk= the Prandtl number 

V \gaLTV) 

In the two-dimensional case, when all variation with respect to y vanishes, these equations 
can be written in component form as (rearranging somewhat) 

aw dP 

at az 

~+eu-vT+2w=<r- 1/2 Rv 2 T 
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d u . 2v dP . 1/2 „ 9 ,„.. 

^- t +eu-Vu-j=-^ + <T l/2 Rv 2 u (16) 

^ t + eU-Vv+~a^Rv 2 v 

dw du_ 
dz dx 

4.2. The Rotation-Dominated System, Q s 

When rotation dominates, the various dimensional scaling factors are related by equating 
the magnitude of the time acceleration term, the Coriolis acceleration and the pressure term in 
eq (9) . Then 

V P 

-=SlV=^ T (17) 

r PqL, 

and these equations imply the following relations between the scaling factors : 

T =l/a, I>= Po LVti. (18) 

Then eq (9) can be written, dropping the prime notation. 

|" + eu-Vu+2i,Xu=-vr+ 2 ^ Ti 2 +EV 2 n 
where 

e=y— =the Rossby number 

E=—j- 2 = the Ekman number. 



(19) 



(20) 



Equation (10) can be written 



dT . _ . 2ATV E 2rr 



(21) 



The magnitudes of the buoyancy term in eq (19) and of the term expressing convective tem- 
perature changes due to the equilibrium temperature distribution in eq (20) arc also equated. 



2gofiJ2ATV 
Vtt~ SQL 

This equation yields the relation between scaling factors 



9 \AT/ 



(22) 



(23) 



Then eqs (19) and (21) can be written 

^+eU-VU+2i 2 Xu=-vP+2^Ti 2 +a 1 'Wu 

~+eU'VT+2Fw=Ea- 1/2 V 2 T. 

where 

F _l /guATV'*_N 
U\ L ) 2U 



(24) 



E=Ea- 1 ' 2 (25) 
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N=l -^-j — ) =the Brunt-Vaisala frequency. 

In the two-dimensional case, when y variations vanish, eqs (24) and the continuity equation 
can be written in component form as 

at ox 

^ t + eU>Vv+2u=v 1/2 Ev 2 v 

^+eU.vw-2FT=-~+<r 1 < 2 Ev 2 w (26) 

ot oz 

^+eu.vT+2Fw=a- 1/2 Ev 2 T 

du dw 

dx dz 

Comparison of eqs (16) and (26) shows them to be identical when the following 
identifications are made: 

System ti s x z u w v T E F 

System Sq z x w u T v R F~ l . 

In general the restriction must also be made that a=l. However, for steady, linear flow, the 
analogy is valid for all Prandtl numbers since it can be incorporated into dependent variables 
by redefinition of these variables. The boundary conditions also must be appropriate for the 
analogy to be complete, and these boundary conditions can be made appropriate. 

In reference 16 several examples of this analogy are presented. One example, waves in 
stratified or rotating fluids, which were discussed in the fully three-dimensional case in section 3 
of the present paper, is shown to be mathematically identical in the two-dimensional case, 
when the identifications above are made. The analogy is much more profound than simply 
one between waves, however, since it applies in dissipative flows with boundaries. Furthermore, 
as noted by eqs (16) and (26), the analogy also applies in rotating flows when stratification 
modifies the flow and in stratified flows with modifying rotational effects. However, as dis- 
cussed by Veronis [16], the analogy breaks down in fully three-dimensional cases. 

5. Nonlinear Considerations 

In this section a very brief account will be given of certain approaches to nonlinear prob- 
lems in stratified and rotating fluids. This is an area of active interest because the equations of 
fluid dynamics are basically nonlinear and because relatively little is known about nonlinear 
phenomena compared with linear phenomena. However, it is also a very difficult area because 
the mathematical tools are limited. 

Only one approach will be examined in any detail. The equations will be derived for finite- 
amplitude, two-dimensional steady flow of a stratified flow over an obstacle. These equations 
are derived and the method of solution briefly discussed because this approach is a unique one. 
Some of the other approaches to nonlinear problems use techniques common to other branches 
of fluid mechanics, and these will only be mentioned. 

The nonlinear equations for steady, two-dimensional flow have been cast into a single 
equation for a stream function by Long [46] and into another equation by Yih [7]. These equa- 
tions can be related simply under certain conditions, but each has been treated in a different 
fashion by its author. The equations will be derived and the method of solution of each briefly 
discussed. It should be noted that both authors have used inverse methods for solving problems. 
Boundar} 7 conditions are selected and a solution is obtained. Then the nonlinear problem which 
was solved is determined as the final step. Such a procedure is not without difficulty and con- 
troversy. There have been several papers in the literature discussing these difficulties [4, 45]. 
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5.1. The Equation of Long for Stratified Flow over an Obstacle 

For a two-dimensional, incompressible flow, the equations are independent of y. Equation 
(2.11) becomes two equations, incompressibility and continuity 

u d P +w d£ =0 (1) 

dx 02 

^+^=0 (2) 

ox dz 

while the equation of motion (2.12) becomes 

/ du . du\ dp /0 . 

p ( u dx +w ^r-dx (3) 

/ dw . dw\ dp , A . 

P K U d- X +W dz) = -dz- dp - (4) 

Here, the motion has been assumed steady as stated above, the equations are written in a 
nonrotating reference frame and the fluid is assumed inviscid so the f=0. 
Equations (3) and (4) can be rewritten as 

d /u 2 +w 2 \ . dp . . n /c v 

d /u 2 +w 2 \ .dp . . n , a . 

p d Z {-^) + I +p(J - pu ^°- (6) 

Here f is the y component of the vorticity, 

du dw 
dz dx 

From these two equations, the pressure is eliminated by cross differentiation; thus 

where q 2 =u 2 +w 2 . 

Equation (2) implies that a stream function can be defined as follows 

- u= Yz w= ¥x (8) 

Then Eq. (2) is identically satisfied. Furthermore, when eqs (8) are used, eq (1) shows that 
the gradient of p is parallel to the gradient of \j/, and therefore 

P= P (+). (9) 

It can be seen that for any functions of these qualities /(p) , g(\f/) 



/ a , a \ [/(") 
\ u a-x +w 



dz 



gin 



df 

dp 

dg 



(4 + »i)£=°- (10) 



In terms of the stream function, the vorticity becomes 

r=-vV. (ii) 
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With these observations, eq (7) can be written 

(«l+»a4)[-|§"+*->>» < 12 > 



noting only that 



dp dp oty dp dp / 9 . d \ 

dx 9 =Ip dx g= T+ Wg= di ( U dx +W dz) gZ 

( d . d\d P 

-\ w Wx +w ¥z)T^ z - 

Equation (12) can be integrated, since the derivative is along a streamline, to give 

-%($+9 z )+P$ = - H *Mf>^ ( 13 ) 

where H*(\p) is a function of \p to be determined from conditions upstream. Substitution for q 2 
and f in terms of \p gives the single, nonlinear equation 

This equation can be written in another form by noting that H* (\p) has a specific form upstream, 
namely 



Thus 



A change of variables from ^ to the upstream variable z , noting that 

d L _dzo_d^_ 1 d 
d\f/ d\// dz U dz 



(16) 



reduces eq (15) to the form 

The transformation (16) can be used only when the upstream velocity U does not change sign. 
In this form, Long [46] observes that the equation becomes linear if pU 2 = constant and 

-j— = constant. 

dz 

Then, with 8=(z — z), eq (17) becomes 

V 2 5+<7 2 <5=0 (18) 



where 

dp 
pU 2 dz 



*>— g - 



Long [46] solved for the flow over obstacles of infinitesimal height first. The resulting 
solution was found to satisfy eq (18) also for obstacles of finite height. However, the streamline 
along the plane z = upstream, which is found to describe an obstacle of finite height, specifies 
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the shape of the obstacle. Subsequent work has extended the analysis to obtain solutions for 
flows over thin obstacles of specified shape and arbitrary height [45]. 

This general method of approach has been extended by Long [45] to other two-dimensional, 
steady, finite-amplitude systems. For example, he has derived a single nonlinear equation for 
the stream function for axisymmetric motion with a body force and for the stream function in 
two-dimensional, steady motion in a rotating fluid. In the latter case, a linear form of the non- 
linear equation can be obtained as described above, and this equation is exactly eq (18) with a 
different expression for a 2 , but with 5 also being the deflection of the stream line from 
its upstream height. 

5.2. The Equation of Yih for Stratified Flow over an Obstacle 

Another single nonlinear equation has been derived by Yih [7] from eqs (l)-(4), and he 
has systematically examined the implications of this equation. In this subsection the equation 
of Yih will be derived. 

Again, eqs (l)-(4) are the starting point for the derivation. Yih [7] observes that the inertia 
effect of a variable density (mentioned in section 2 of this report) can be accounted for, defining 
new velocity variables n' , w' , as follows: 



(19) 



V Po 

w' = -\ — w 

Y po 

where p is a constant reference density. Then Eqs (l)-(4) become 

/ , du' , du'\ dp /ooN 

/ ,dw' , dw'\ dp , 00 , 

(Note that these equations follow immediately from the observation made in eq (10) with 

/(p)=Vp/poO 
Equations (22) and (23) can be rewritten as 

ozj 2 ox 



Po( lif Q x + W 



( , d . , d\w' 2 , dp , , OK v 

p { u w x +w d-zh=~ w d- z - gpw (25) 

when eq (22) is multiplied by u' and (23) by w' . Addition of eqs (24) and (25), with the observa- 
tion that 



yields the relation 



„»'-»(«' £+»' l) z= ("' 1+"' m) 9 "' 
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This relation implies that the quantity 

H= Po - Y-p+gpz, 

known alternately as the Bernoulli quantity [7] and the "total head" [4], is a constant along a 
stream line. It will, however, vary from stream line to stream line. Yih [7] defines a stream func- 
tion ty 1 based upon the velocities u' and w f as follows: 

Then eq (21) is identically satisfied, and since H is a function depending upon the stream line, 

With the definitions and observations made above, the single nonlinear equation of Yih 
can be derived. Equations (22) and (23) can be written as: 



d /u' 2 +w' : 



)+|j+p oW T=0 (28) 



Po 3* V 2 

Po|(^>g + P,-Po.Y=0 (29) 



where 



,_6V 6V 
"" dz dx 

is the vorticity based upon the velocities u', w'. Multiplying eq (28) by dx and eq (29) by dz 
and adding yields 

r 7/ '2_i_ 7/ ,/2 -i 

d [po g +P+pgzj-Po{'dt'-zgd P =0 (30) 



where 



Since 



d=dx^ — Ydz-Tr 



ff=Zr(^ , )andf , =vV / , 



eq (30) becomes 



0tlx zqdp 1 dH _ /01 v 

VV +— 377 ^ =0 ^ ( 31 ) 

Po <# Po a^ 



and this is the equation of Yih. 

Yih systematically examined solutions to this equation when 

dp i dH 

d\p' d\f/' 

are linear functions of \[/' ; then eq (31) is linear and classes of solutions to the original nonlinear 
equation are found. For details of these solutions, see references [7]. 

5.3. Other Nonlinear Problems 

The purpose of this subsection will be to mention some other nonlinear problems in strati- 
fied or rotating fluids. In some cases the mathematical development is exactly the same as in 
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other areas of fluid dynamics; in others the mathematical techniques have been developed for 
problems in different areas and have subsequently been applied to flows of stratified or rotating 
fluids. 

The free surface flow of a shallow layer of homogeneous fluid over an obstacle is an example 
of a stratified (layered) flow in which the governing equations are exactly the same as the one- 
dimensional, compressible-flow equations for a poly tropic gas in which the gas pressure is 
proportional to the density squared. This analogy has been known for some time [47], and the 
governing equations are hyperbolic [26] as defined in section 3 of this report. Perhaps the most 
interesting aspect of this analogy is that an identification can be made between the internal 
Froude number in the stratified flow and the Mach number in the compressible flow [7, 45, 47]. 

The internal Froude number is defined in eq (2.5) ; here it is rewritten slightly as 

F= u 

V/ 2 (32) 



(*W 



where U is the uniform upstream velocity in the shallow layer of fluid, p { is the density and L 
the thickness of the layer of fluid and p 2 is the density of the fluid overlaying this layer. The 
velocity 



(^r 



is that of infinitesimal long waves on the layer of fluid. The Mach number is defined in eq 

(2.10) as 

M=U/c 

where U is the characteristic flow velocity and c is the speed of sound waves. 

When the Froude number is loss than one, the flow is said to be subcritical; when /'\>1, 
supercritical and when F=l, critical. In each of these cases, the surface elevation of the flow over 
the obstacle changes in character. In subcritical flows, the surface elevation decreases from its 
upstream value as the fluid flows over the obstacle. In the supercritical case, the surface ele- 
vation increases as the flow passes over the obstacle. When the flow becomes critical, phenomena 
known as hydraulic jumps or bores occur. For a description of the flow regimes as a function of 
internal Froude number, see references [45] and [4]. 

The point to note here is the analogy between stratified flow and compressible flow. The 
Froude number is the parameter which delineates various regimes exactly as the Mach number 
delineates regimes of flow in gasdynamics. The phenomenon of the hydraulic jump is exactly 
analogous to the shock wave in gasdynamics. Furthermore, for a variety of flow conditions, the 
mathematical formulation and solutions are identical. 

An area of rather intense study presently is that of nonlinear waves [26, 48, 49] ; nonlinear 
wave phenomena appear in many fields, and much of the mathematical apparatus is applicable 
to all of these fields. Formation of shock waves in compressible flows or hydraulic jumps in 
water-wave theory is an example of a nonlinear wave arising from a hyperbolic equation. More 
recent attention has been focussed on nonlinear, dispersive waves, where the nonlinear tendency 
to steepen is balanced by the dispersive tendency to spread them. Examination has been made 
of both periodic (cnoidal) waves and nonperiodic (solitary) waves in such systems. 

Another particularly active problem area in the theory of nonlinear, dispersive waves is 
that of wave-wave interactions. Conditions are examined for which the nonlinear terms in the 
equations of motion can permit resonant coupling between wave modes of the linear system. 
One example where such coupling can occur is in the interaction between internal and surface 
waves [49]. Another example, which has been intensively investigated recently, is that of 
coupling between laser radiation and oscillation modes in a plasma; these so-called parametric 
instabilities have been examined for possible application in laser-driven fusion [50]. 

As final examples of nonlinear problems, boundary layers in rotating fluids [9] and non- 
linear convection in stratified fluids [4] are mentioned. In the next section of this report some 
of the linear features of these problems will be discussed. 

6. Transport and Diffusion Processes 

In this section some interesting problems will be discussed related to transport and diffusion 
processes in a stratified or rotating fluid. Since the work in this area is immense, only three 
problems have been selected for discussion, and even these are discussed only briefly. However, 
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it is felt that this presentation will convey the type of studies being carried out. Reference 
should be made to the monographs and papers cited in this section. 

The three problems to be discussed are contained, rotating fluids, convection in stratified 
fluids, and turbulence in stratified or rotating fluids. The problems of contained, rotating fluids 
has been reviewed extensively b}^ Greenspan [9], who has also been a major contributor to 
the research in this area. A more recent review of the problem of "spin-up," which is a result 
of viscous effects and boundaries, has been given by Benton and Clark [14]. The analogy 
between rotating and stratified fluids [16], discussed in section 4, allows some of the results 
from rotating fluids to be carried over to stratified fluids. Convection in stratified fluids covers 
a large number of phenomena. Turner [4], who has been a major contributor in this area, has 
recently reviewed these phenomena. Finally, turbulence, which in many respects acts as a 
diffusive effect, in stratified fluids will be discussed. Turner [4] has recently described some 
important features of turbulence in stratified fluids. Turbulence in the atmosphere has been 
discussed by Lumley and Panofsky [55], and Phillips [5] has reviewed turbulence in the ocean. 
Work in this very important area has been growing rapidly. 

6.1. Contained Rotating Fluids 

The most comprehensive review of flow phenomena in rotating fluids is the monograph 
of Greenspan [9]. In that book he presents a systematic development of the formulation and 
solution of problems concerned with contained rotating fluids. Boundaries and viscous boundary 
layers produce such a dominant effect in these fluids, they are discussed briefly here. First 
the basic linear equations are derived. Then the Ekman layer, a fundamental boundary layer 
in rotating fluids, is discussed. Finally, the phenomenon known as spin-up is mentioned. 

The basic equations for the analysis of contained rotating flows are eqs (2.11) and (2.12). 
The fluid is assumed to be homogeneous, incompressible and viscous; therefore, eq (2.11) 
becomes p= p and 

V-u=0. (1) 

Equation (2.12) becomes 

Po ^+2p OXu+vP=Po^V 2 u. (2) 

Here P is the reduced pressure defined earlier, the body force due to gravity has been dropped 
because the fluid is homogeneous and the force f in eq (2.12) has been taken to be due to 
viscosity. 

The characteristic length, time and velocity scales, L, ft -1 and U are chosen to nondimen- 
sionalize the variables. Therefore r, t, u, Q, = tti z and P are all scaled by these quantities to 
become, Lr, Q~H, Uu, £2i 2 and p ^ULP, and eqs (1) and (2) become in dimensionless form 



V-u=0 (3) 

du 
where 



n +€U.Vii+2i 2 Xu+vP=#V 2 u (4) 



is the Rossby number and 



U_ 
6 QL 



E= ^U 



is the Ekman number. When the equations are linearized, €=0, and eq (4) becomes 

^+2i 2 Xu+vP=Ev 2 u- (5) 

Equations (3) and (5) can be used to derive a single equation for the pressure. When the 
curl of eq (5) is taken, the equation for the vorticity co=vXu is obtained 

(|-^)vxu= 2 |a. (6) 
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When the divergence of eq (5) is taken, again using eq (3) , 

V 2 P=2i 2 .(vXu). (7) 

The operator 

<£-*■) 

is applied to eq (6) and eq (5) is used to obtain 

(|-^V 2 ) 2 (VXu)=2^[-vP-2i z Xu] (8) 

Then the z component of this equation, using eq (7), gives the desired equation for the pressure. 

/ d V d 2 P 

(ii- E V+ p --*w- (9) 

Some idea of the difficulty of problems of this type is obtained by examination of eq (9). 

A very fundamental role in the theory of rotating flows is played by the Ekman layer. 
In this layer a balance is struck between viscous and Coriolis forces. In the simplest case, this 
balance can be seen by the following model. Consider a half-space of fluid 2>0, rotating with 
constant angular velocity ft. Assume that a uniform velocity u=U ) v=0 is imposed at the 
surface 2 = 0, and look for a steady-state solution which is independent of the horizontal co- 
ordinates x and y. The scaled eqs (5) become 

~ 2v = E d?' 2u = E d?- (10) 

The solution to eqs (10) which satisfies the boundary conditions is 

M=exp (~A) C0S A 

(id 



;=exp (-jb) sin 



z 



From the continuity equation (3) and the third of eqs (5) it is found that the vertical velocity 
w=0 and that the pressure is constant. 

Two important features of the Ekman layer can be seen from the solution (11). First, the 

Ekman layer is a boundary layer of thickness -*/— ; the velocity, which is unity (in nondimen- 

sional units) at 2=0, decays exponentially over this thickness. Second, the velocity vector 
rotates as z increases. Plotted in polar coordinates, the velocity forms a spiral shape, the so- 
called Ekman spiral [9, 16]. 

Examination of the dynamics of the Ekman layer under more realistic conditions shows 
that the vertical velocity is not zero as determined by the simple model above [9]. Rather, the 
vertical velocity increases exponentially through the Ekman layer from zero at 2=0 to a 
value of — E 1/2 . When Sis small, as is usually the case, this vertical velocity is small. However, 
this suction velocity dramatically affects the inviscid flow. 

Spin-up is one problem in which the Ekman layer plays a fundamental role. Consider a 
container filled with a homogeneous fluid and rotating with the fluid about a vertical axis at 
a constant angular velocity. Then the angular velocity of the container is suddenly increased 
to a new value and maintained at that value. Along the top and bottom surfaces of the con- 
tainer Ekman layers form. In each layer the centrifugal force is too great for the prevailing 
pressure gradient, and the fluid is driven radially outward. By continuity, the suction velocity 
into each layer arises, and this velocity produces secondary circulations within the contained 
fluid, including additional boundary layers of a different nature along the vertical walls. Ulti- 
mately (after a time scale L/v 1/2 tt 1/2 , determined by the length scale L and the Ekman suction 
velocity E l/2 QL) these secondary flows lead to a state of solid body rotation of the fluid at the 
new angular velocity. This rather complex process is called spin-up [9, 14]. 
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6.2. Convection in Stratified Fluids 

In a recent book, Turner [4] has given a review of many interesting and practical con- 
vective phenomena in stratified fluids. In this subsection, following Turner, a very brief dis- 
cussion of these phenomena will be presented. It will address convection from isolated sources 
of buoyancy, convection from heated surfaces and so-called double diffusive phenomena. 

a. Buoyant Thermals and Plumes 

In many situations the buoyancy in a fluid can supply the energy for the mean motion 
within the system. An example of such a situation is the thermal, well known to meteorologists 
and pilots, in which a mass of warmer, lighter air rises due to buoyancy in cooler, more dense 
surroundings. Another example is the thermal plume created when hot water from a power 
plant is discharged into a waterbody, the discharge water rising in the cooler ambient water. 
(Effluents from chimneys in the atmosphere also often fall into this category when the effluent 
is more buoyant than the surrounding air.) 

These two examples display the properties of two broad categories into which flows of this 
nature have been divided by Turner [4], The first, called the thermal, is characterized by a 
sudden release of buoyancy within a fluid. The buoyancy remains confined to a limited amount 
of fluid, which, as it rises, loses connection with the source that produced it. The second, called 
a plume or buoyant jet, is characterized by a continuous source of buoyancy. Rising buoyant 
fluid is always connected by additional buoyant fluid to its source. 

The driving force in both of these types of motion is the buoyant force. Lighter fluid will 
be buoyed upward in its surrounding heavier fluid. (Likewise, heavier fluid will fall in lighter sur- 
roundings.) This fluid will rise causing a mean convective motion and smaller scale motions 
called entrainments near the interface between the buoyant fluid and its surroundings. Gen- 
erally, the entrainment is a turbulent process in which surrounding fluid is captured into the 
thermal or plume fluid. Entrainment causes the volume of buoyant fluid to be increased, but 
the net buoyancy of each fluid element to decrease relative to the surroundings. 

The mathematical techniques, in addition to numerical ones, generally used to analyze 
buoyant thermals and plumes are similarity methods. A major consideration in the analysis 
of these flows is the entrainment and mixing of the ambient fluid, since entrainment determines 
both the mean size and the net upward flow velocity of the thermal or plume. Other complicat- 
ing features in the analysis are nonhomogeneous ambient media and an ambient steady or 
unsteady mean cross flow. The steady cross flow causes a plume to be "bent over" in the direc- 
tion of flow, changing the mean properties of the plume and often the character of the entrain- 
ment as well. 

b. Convection from Heated Surfaces 

A fundamental problem with a long history in the theory of convective flows is the Benard 
problem [52]. Fluid, contained between two horizontal surfaces, is heated from below. The 
increased buoyancy tends to destabilize the fluid, and when the dissipative effects of thermal 
conduction and viscosity are overcome, the fluid will overturn. The parameter which governs 
this stability problem is the Rayleigh number [4, 52] 

Ra=S»§^. (21) 

Here, all of the symbols have their usual meanings. AT is the temperature difference across the 
fluid layer of thickness d and K is the thermometric conductivity. The results of the experi- 
ments and analysis of this problem show two main features. First, a critical Rayleigh number 
must be exceeded before instability occurs; this implies that a certain temperature difference 
must be exceeded for a given fluid and layer thickness. Second, when convection sets in, the 
motion is organized in regular, steady cellular patterns. Chandrasekhar [52] has given an 
encyclopedic review of the linear aspects of this instability problem, including the effects of 
rotation and magnetic fields. Subsequent work has concentrated on finite-amplitude aspects 
of the convection [4]. 

Turner [4] also reviews other studies of convection caused by heated surfaces. These 
include intermittent thermal or plume formation at a horizontal surface, so-called penetrative 
convection (where convection initiated in an unstable region penetrates into an adjacent stable 
region), and convection layers induced at vertical or sloping boundaires. 
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c. Double-Diffusive Convection 

An area of study in convective phenomena, which is of relatively recent origin, is double- 
diffusive phenomena [4]. A fluid is considered for which the density of the fluid is determined 
by the concentration of two properties, which have differing diffusivities. In the ocean, for 
example, the temperature and the salinity are the two properties; each affects the density of 
the water by its concentration, and the salt and temperature have different diffusivities. 
These phenomena appear to be applicable to more problems than only oceanic ones. 

In the ocean, the concentration of salt and the temperature of the water determine its 
density. An increased salinity increases the density whereas an increased temperature de- 
creases it. Therefore, the possibility of interesting double-diffusive instabilities arises. For 
example, if colder, less saline water underlies warmer, more saline water, such that the density 
of the lower layer is greater, then the fluid should be stable. However, if a fluid parcel from the 
lower layer is displaced upward, due to the difference in diffusivities, the parcel heats up while 
retaining its lower salinity. Then the parcel is less dense than the surrounding fluid and rises, 
leading to an unstable situation. This instability is the mechanism for a "salt fountain, " which 
is the flow through a heat-conducting pipe of the colder, less-saline fluid to the top of the over- 
lying layer of warmer, saline water (or the opposite flow) . 

When cold, fresh water is above hot, saline water, such that the density increases down- 
ward, another type of instability can occur. A displaced fluid element of hot saline fluid will be 
more dense than its surroundings and will tend to sink. In the absence of any diffusion processes, 
the fluid element would oscillate at the local Brunt-Vaisala frequency. However, because of 
the relatively large thermal diffusion, the displaced element would cool somewhat before 
reaching its original equilibrium position, and, therefore, would plunge lower than without 
diffusion. Heating during the lower portion of the oscillation cycle would then cause the fluid 
parcel to rise higher during its upward motion. Repeated cooling and heating during the 
oscillation cycle causes oscillations of growing amplitude, a process known as overstability. 

These two instabilities can lead to interesting phenomena known as salt fingers and micro- 
structure. Salt fingers are long vertical columns of warm salt water in cooler, less saline water; 
they are believed to result from the former instability. Microstructure is a series of horizontal 
layers of fluid of relatively uniform properties separated by rather abrupt changes in these 
properties. It is believed that some microstructure is a result of double diffusion convection [4]. 

6.3. Turbulence 

Turbulence is one of the most prevalent phenomena within fluids. Its motions are often 
mainly responsible for the transport and diffusion of mass, momentum and energy within these 
fluids. Because of this prominent position, it has been examined extensively. But, because of its 
inherently nonlinear and random nature, progress in understanding and modeling it has been 
slow. In this subsection no attempt will be made to review the extensive literature in this area; 
rather, a very brief discussion will be given to indicate its importance in the study of stratified 
and rotating fluids. 

Turbulence in stratified fluids is an especially interesting and challenging problem. Mixing 
of a stably stratified fluid requires the expenditure of energy, and, therefore, vertical turbulent 
mixing is inhibited by stratification. The parameter which governs this interaction is the local 
Richardson number. The local Richardson number compares the suppression of turbulence by 
the mean buoyancy forces to the production of turbulence by Reynolds stresses [5]. 

Turbulence in the atmosphere, which is generally stratified, has been and continues to be 
an area of active interest [53-55]. The atmospheric boundary layer is also influenced by rota- 
tional effects [56]. In fact, a model of this boundary layer is an Ekman layer in which Coriolis 
forces are balanced by Reynolds stresses (rather than viscous stresses as described earlier). 

Similarly, turbulence in the ocean and in lakes is a dominant influence on their struc- 
tures [5]. The interaction between turbulence and buoyancy gradients is often manifested by a 
rather sharp interface between a turbulently mixed region and an adjacent stably stratified 
region [4, 5]. It is this interaction which has been used to explain the distinctive temperature 
structure in temperate lakes, for example [57]. 

For many purposes in geophysical applications, including the two specific ones cited 
above, turbulence acts as a diffusive process. Then the effects of turbulence are modeled by 
mixing-length theories, where the nonlinear, random mixing is related to the mean-flow gra- 
dients through turbulent diffusion coefficients. For discussion of the modeling and the effects 
of turbulence, the references cited in this section should be consulted. 
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7. Some Applications 

The basic concepts in studies of stratified and rotating fluids have wide applicability, 
particularly in the area of geophysics. Some of the applications have already been mentioned 
in the introduction and at other parts of this report. In this section, a brief summary of the 
geophysical applications is presented. This summary will not be complete, but will convey the 
rich variety of interesting problems and results. 

7.1. Physical Oceanography 

In the ocean a very large variety of fluid motions occur, ranging in length scale from milli- 
meters to tens of thousands of kilometers. It is the largest scale circulations which differ most 
from classical fluid flows because curvature and rotation of the earth, stratification produced 
by solar heating and even the differential gravitational effects producing tides dominate these 
motions [58]. However, even the smaller scale motions are profoundly influenced by rotation 
and stratification. 

The ocean is a very thin layer on a nearly spherical globe. In the deepest region, the 
Marianas Trench, the ocean is about 12 km deep. The radius of the earth is about 6400 km, 
so that the ocean depth is less than 0.2 percent of the radius of the earth [59]. For consideration 
of many fluid motions the curvature of the earth can be neglected; but for larger curculations, 
the curvature must be considered. A common approximation used to simplify anatysis of 
oceanographic circulations is the "/3-plane approximation. " The fact is used that the depth A 
of the ocean is much less than the radius R of the earth. It is assumed that the lateral length 
scale L for motions is much greater than the depth and that Z<<i?, the most questionable 
part of the approximation [58]. Then the equations of motion for the sheet of fluid on the rotating 
globe are reduced to equations for a thin sheet of fluid in a rotating tangent plane. 

The waters at the equator are heated by solar radiation and those at the poles undergo 
net radiational cooling [59]. In the absence of rotation, therefore, the more buoyant surface 
fluid near the equator would spread toward the poles. However, the Coriolis force curves this 
spreading, producing large clockwise eddies in the northern hemisphere and counterclockwise 
eddies in the southern hemisphere. These large circulations are also known as g} r res. 

Both heating and salinity determine the density and therefore the buoyancy of the fluid. 
The temperature and salinity combine generally to form a very distinctive density profile in 
which an approximately uniform-density layer overla}^s another layer of slowly varying prop- 
erties. The narrow region which separates these layers, where the temperature varies rapidl} r 
and the Brunt-Vaisala frequency reaches a maximum with depth, is called the thermocline. 
This stratification has important consequences for oceanic circulations. 

Another very important influence on the motions of all scales in the ocean are the winds 
[60, 9]. The large scale, nearly steady atmospheric circulations are discussed in the subsection 
on meteorology. They are driven by solar heating and the earth's rotation in much the same 
way as the oceans. These winds transfer mechanical energy to the ocean surface by generating 
currents and waves and receive heat, primarily through evaporation from the ocean. The 
large-scale oceanic circulations are driven in a major way by this interaction [60]. 

The land masses also have a major influence on the oceanic gyres. The continents act as 
boundaries inhibiting the free circulations. The boundary layers at these surfaces influence 
the overall fluid motions in important ways [61]. On smaller scales, the shape of land masses 
cause many interesting phenomena such as wave diffraction and wave breaking. 

The gravitational influence of the sun and moon also affect the largest scale motions in 
the ocean, forming the well-known oceanic tides. Tidal theory considers long wavelength 
waves in shallow water on a rotating globe and has been reviewed recently by Hendershott 
and Munk [62]. 

In addition to these large scale motions, many smaller scale motions occur. Both internal 
and inertial waves occur within the ocean [5, 6] and surface waves are created by the wind 
[5, 47]. Also, at the surface, evaporation and condensation are important interchanges with 
the atmosphere, altering the local salinity and providing a major source for energy transfer 
with the atmosphere [60]. 

Finally, the pervasive influence of turbulence is felt. In the upper ocean above the ther- 
mocline, it is a major cause of mixing and diffusion [5, 60]. Likewise, near all land masses, 
turbulent boundary layers are dominant [61]. Most of the upwelling motions, the rise of nutrient- 
rich lower waters, are influenced by turbulence. Even the deep oceanic currents are subject 
to its nonlinear, random processes. 
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7.2. Physical Limnology 

Many of the dynamical considerations associated with oceanography are also associated 
with limnology. The differences arise primarily in the large scale motions, where lake circula- 
tions are determined by the confines of the basin. Lakes are also generally composed of fresh 
water, but even this is not always the case. 

Compared with the differences, the number of dynamical similarities is great [19]. Rota- 
tion affects the dynamical behavior of larger lakes [19, 13]. As in the ocean and the atmosphere 
the horizontal dimensions of motions are much larger than the vertical dimension. Lakes are 
bounded by the atmosphere and derive mechanical energy for driving circulations from wind 
frictional effects. Thermal energy is provided by the sun, and thermal interchange takes place 
with the atmosphere. Finally, as in the ocean and the atmosphere, turbulence is a dominant 
influence on the motions. 

In the remainder of this subsection some of the interesting features of modeling lake 
motions will be discussed. As with the other subsections on applications, no attempt can be 
made to be complete. 

The Reynolds number determines when a flow is laminar and when it makes a transition 
to turbulence : for a Re3molds number greater than a critical value, the flow is turbulent. For 
a length scale equal to the depth, the velocity must exceed a minimum value for turbulence. 
For a very shallow lake of 100 cm, the flow velocity must exceed 0.03 cm/s [19]. Clearly then, 
the turbulent state is the natural one in a lake. 

Because turbulence is a random process, statistical averages in time or space are taken to 
describe the character of the turbulence. This averaging brings up a very fundamental question 
of separating turbulent motions from longer-time-scale (or length-scale) motions associated 
with mean flows [55]. Assuming that this separation can be made, the nonlinear nature of the 
turbulence brings up a second fundamental problem, the so-called closure problem. If statistical 
averages of the equations of motion are taken, the nonlinear terms always couple lower order 
statistical correlations to higher order ones. 

In applications such as limnolog t y (and oceanography and meteorology), it is generally 
assumed that second order correlations are related to gradients of the mean flow through 
coefficients of eddy diffusivity [19]. When this is done, the effects of turbulence arise in the 
equations of motion as additional diffusion terms. A major problem, then, is to determine both 
the size and the variation of these coefficients with position and properties of the flow [19, 57]. 

The primary source of mechanical energy for currents in lakes is the wind, introducing 
both mean circulations and mixing into the main body of the water. The wind-driven circula- 
tions begin with water movement in the downwind direction, but since the lake is enclosed, a 
return current is also induced. Both influents and effluents can also affect circulation patterns, 
however. These circulation patterns are also profoundly affected by rotational and stratification 
effects. 

In larger lakes the equations of motion and many of the approximations used in analyzing 
these equations are the same as those used in oceanography. In particular, geostrophy and 
Ekman layers are very important in interpreting limnological flows [19, 13]. 

Also, stratification is a dominant influence in the hydrodynamics of lakes. The stratification 
in lakes arises primarity from surface heating, and, therefore, lakes undergo an annual cyycle in 
overall stratification [19, 57]. During the spring and summer when the solar and atmospheric 
heat input is largest the degree of stratification, or the stability of the lake increases, and 
during fall and winter the lake cools and destabilizes. 

Stratification affects both the current in the lake and the waves which the lake can sustain. 
In the absence of stratification currents in a lake would vary smoothly from a maximum value 
near the surface to zero at the bottom. In the presence of stratification the flow velocity will 
generally reverse in direction as a function of depth. Such behavior is very important in ana- 
lyzing the mechanics of lakes. 

As in the ocean a thermocline is often formed within a lake. Because of the seasonal varia- 
tion in the stratification, the depth and thickness of the thermocline vary considerably. Because 
of the thermocline, however, a lake is often modelled as two superposed layers of fluid of 
different density, the so-called two-layer model. Internal waves, as well as surface waves, can 
be sustained within a lake due to the stratification. 

Because a lake is confined by its basin, it has normal modes of oscillation. The term seiche 
is used to describe stationary oscillations of a lake or a large independent portion of a lake. 
Seiches have been studied extensively. Both stratification and rotation have important effects 
on seiches [19, 13]. 
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7.3. Meteorology 

The atmosphere also is a very thin layer of fluid on the almost spherical earth, the thickness 
of the significant portion of it being about 10 km, or a fraction of 1 percent of the radius of the 
earth. Therefore, horizontal dimensions of motions are generally much larger than the vertical 
dimension. The stratification of the atmosphere enhances this behavior since it tends to inhibit 
vertical motions. 

In an interesting review of meteorology, Sutton [63] concentrates on three systems of 
motion which are determined by " their scales of horizontal length and their mean 'lives'. " The 
largest motions are "huge quasi-permanent wind belts of the general circulation." The inter- 
mediate ones are the "smaller, short-lived cellular-like systems (highs and lows) that determine 
the day-by-day occurrence of weather. " The smallest and shortest-lived motions are "the 
transient eddies that are responsible for turbulent transfer processes in the earth's boundary 
layer." 

The general circulation is determined by the same large-scale effects that determine 
oceanic circulations, namely rotation, curvature of the earth, heating at the equator and cooling 
at the poles. These large scale circulations are associated with climate, the monthly or seasonal 
properties of the atmosphere averaged over long periods of time. Stewart [60] describes the 
atmosphere as a heat engine, receiving solar heat from the surface of the earth and the ocean, 
converting a portion of it to mechanical energy and discharging the remaining heat at lower 
temperatures. As mentioned in subsection 7.1 the wind of the atmosphere drives the oceanic 
circulations. In turn the ocean supplies a major fraction of the heat which the atmosphere 
uses in producing mechanical energy, some through direct heat transfer to the atmosphere, 
but most through evaporation at the surface of the ocean. The evaporation-condensation 
cycle is a major source of energy transfer between the ocean and atmosphere, which are very 
intimately linked dynamically [60]. 

Models for weather prediction have been reviewed by Sutton [63] and by Phillips [64] and 
include the effects of the earth's rotation, curvature of the earth and stratification of the 
atmosphere. Often the vertical equation of motion is replaced by the equation for hydrostatics, 
the so-called hydrostatic system. A further approximation sometimes made is that the horizontal 
equations can be replaced by a balance between pressure and the Coriolis forces, the geostrophic 
relation. An energy equation to account for the important energy fluxes in the atmosphere 
is included, and account must be taken of the water-vapor content of the air. The full equations 
or the various approximate sets are then integrated numerically, often over a 24-hour period. 
Extended forecasts have shown that a small predictive capability exists beyond 10 days [64]. 

Momentum and energy transfer between the earth and the atmosphere are important 
in about the lower 1000 meters of the atmosphere, the atmospheric boundary layer [53-56]. 
In this region the turbulent diffusion is balanced by the Coriolis force so that this layer is an 
Ekman layer. Radiative energy absorption takes place primarily at the surface of the earth 
so that atmospheric heating occurs from the bottom. (This is to be contrasted with the ocean 
where heating takes place at the top.) Therefore, thermal energy transfer and buoyancy effects 
are dominant in the atmospheric boundary layer. 

The study of the turbulent boundary layer is called micrometeorology. Important fea- 
tures of this study are the small scale, turbulent mass, momentum and energy transfer. Micro- 
meteorology has become a prominent area of study as concern has increased about diffusion 
of pollutants in the lower atmosphere. 

7.4. Thermal Discharges 

A major problem of practical concern in recent years is that of discharges into the atmos- 
phere or into bodies of water. In subsection 6.2a on Buoyant Thermals and Plumes, some 
fundamental fluid-mechanical features of buoyant discharges (usually produced by heating 
the effluent and therefore called thermal discharges) and techniques for calculating these 
features were reviewed. Here, some general comments on the behavior of thermal discharges 
and the practical considerations related to these discharges will be made. 

The mixing between a thermal discharge or plume and the ambient fluid can be viewed 
to take place in two parts. Near the outfall, the entrainment and mixing is dominated by the 
turbulence induced by the discharge, and this region is called the near field. In the far field, 
the entrainment and mixing are dominated by the naturally occurring turbulence in the ambient 
fluid. 

The near field has been discussed in a variety of papers [4, 65-68]. If the effluent is dis- 
charged with a substantial initial momentum, the effects of buoyancy can often be neglected 
in the immediate vicinity of the outfall; then standard results on the growth rate, entrainment 
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and mixing can be used from the theory of turbulent jets. When the ambient fluid is flowing 
at right angles to the discharge, the plume will quickly be bent over in the direction of the 
cross flow. However, the initial buoyancy and momentum will continue to be important in 
the dynamics of the plume. 

In the far field, diffusion is accomplished by the ambient turbulence, although the initial 
momentum and buoyancy may still be important. When the ambient fluid is stably stratified, 
the diluted effluent will reach a level where its density equals that of the ambient fluid, and 
will spread almost entirely in the horizontal direction. The spread of a passive material in a 
flowing or stratified, turbulent medium has been reviewed extensively [69-71, 5]. 

The work on plumes is divided into two groups by application, atmospheric plumes and 
discharges in water. Each application has particular, interesting problems and an extensive 
literature associated with it. For example, in the atmosphere the initial buoyancy is deter- 
mined not only by the temperature, but also by the average molecular weight and the particu- 
late content of the discharged gas. Whereas the gas will generally rise and mix as described 
above, the particles will settle out at a rate dependent upon their mass and size distribution. 
When water vapor is present, condensation, with an associated heat release and increase in 
buoyancy, can take place. In addition, chemical reactions, and particularly photochemistry 
(induced by the sun), arc important in consideration of atmospheric pollutants. Therefore, 
the equations describing the chemical kinetics must be carried along with the equations for 
the plume dynamics. Finally, in engineering design, stack height and placement are important 
for consideration of the dispersion of effluent material. 

Similar problems must be considered when thermal discharges in water are analyzed. 
However, some of the pioblems and the emphasis of the problems are different. In a water 
body, a surface is piesent and this can alter the dispersion dynamics greatly. In addition, as 
described earlier, a distinctive thermal structure including a thermocline, is often present, 
trapping the discharge at a particular level and causing two-dimensional, horizontal spreading. 
If the discharge does reach the surface, the spreading will be two-dimensional and heat will 
be lost to the atmosphere directly and through evaporation. From the engineering-design 
viewpoint, the position of the water intake and discharge and the design of the diffuser at the 
discharge are important. 

8. Some Important Problem Areas 

In this section some important problems related to stratified and rotating fluids will be 
discussed. The problems can be broadly classed into four categories, turbulence, nonlinear 
processes, numerical computation of flows and applications. Each category will be discussed 
briefly; as in preceding sections, no attempt will be made to be complete. The emphasis will 
be on theory, but it should be noted that proper understanding arises only when theory and 
experiment supplement each other. 

Turbulence remains one of the most difficult and challenging phenomena in fluid mechanics. 
Because it is a nonlinear, random process, the progress in understanding and analyzing it has 
been slow. Some years back there was a great deal of research activity and several books 
summarizing the status. Some of this work concerned homogeneous, isotropic tuibulence and 
some turbulent shear flows. 

Recently, there has been an increase in interest again, and the emphasis has now been 
on turbulent shear flows and on turbulent mixing in the presence of stratification. Both of 
these aspects of turbulence enjoy wide applicability, as seen in this review. They also have 
another feature in common: the mean fluid properties are anisotropic. In shear flow, a direction 
is determined by the vorticity of the mean flow. In stratified fluids, a direction is determined 
by gravity. The anisotropy might be expected to affect the turbulent eddies in fundamental 
ways: eddies with their vorticity parallel to the direction determined by the mean properties 
will be influenced differently than those perpendicular to the mean properties. With the re- 
surgence of interost in turbulence, it is anticipated that many of the problems associated with 
turbulence in an anisotropic medium will be addressed, and progress in understanding it will 
be made. 

Early in this report it was noted that linear problems have received much more attention 
than nonlineai ones because the tools of linear analysis are more highly developed. The analysis 
of nonlinear problems has always required elements of ingenuity, patience and luck, and in- 
dividual problems generally require unique approaches. However, an increased number of 
solutions to nonlinear problems has increased the variety of available approaches to these 
problems. Some approaches can be used directly or generalized to attack other problems. So, 
many areas of nonlinear analysis are developing. In particular, nonlinear waves, convection 
and motions in rotating or stratified fluids are all being examined to an increasing extent. 
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Continued effort in these areas is needed. In addition, the methods of analysis should be classified 
and generalized to aid the attack on other nonlinear problems. 

With the advent of high-speed computers, numerical solutions to a wide variety of physical 
problems have become a reality. As a result, computational physics and computational fluid 
mechanics have become established branches of theoretical physics and theoretical mechanics. 
Finite difference techniques have been developed to integrate ordinary differential equations 
and initial-value and boundary-value problems associated with partial differential equations. 
Although enormous effort has been expended in recent years in developing and increasing the 
efficiency of numerical techniques, more is needed. With the guidance of caieful analysis, the 
computation of fluid flows, especially stratified and rotating flows, should be pursued vigorously. 

Finally, the application to physical problems should be made of the basic results deter- 
mined for stratified and rotating fluids. Much work has already been performed along this 
line, and some of it was discussed in the preceding section. However, many problems still 
remain in the earth sciences, in astropr^sics and in more technologically oriented areas such as 
those discussed in section 7.4. 

9. Summary 

A survey has been presented of phenomena in stratified and rotating fluids. Care was 
taken to define important quantities, to discuss basic concepts, to derive the fundamental 
equations and to present the basic nondimensional parameters associated with these flows. 
After a rather extensive section on waves, the analogy between stratified and rotating flows 
was discussed. Then nonlinear processes, and transport and diffusion processes were reviewed. 
Although this report is rather brief in parts, it displays the rich variety of phenomena in strat- 
ified and rotating fluids. It also tabulates many of the important references. 

This report also contains a discussion of some applications. Physical oceanography, physical 
limnology and meteorology are all areas of application in the earth sciences. The variety of 
applications was stressed at the expense of depth and completeness. The more technological 
application of plumes in bodies of water or in the atmosphere was also discussed. 

Finally, some important problem areas were briefly surveyed. These areas are turbulence, 
nonlinear processes, numerical computation of flows and further applications. 
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